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INTRODUCTION 


Nuclear  Astrophyslcal  Interest 


TO  1  fc\ 

The   C(a,y)   0  reaction  is  of  considerable  interest 
in  nuclear  astrophysics,  especially  in  studies  of  stellar 
nucleosynthesis  and  models  of  stellar  Interiors.   This 
reaction  occurs  in  the  helium  burning  sequence  and  has 
a  direct  influence  on  the  relative  abundances  of  oxygen 
and  carbon.   The  helium  is  produced  from  hydrogen  burning 
by  the  proton-proton  and  carbon-nitrogen-oxygen  cycles  (1,  2) 
After  the  hydrogen  Is  consumed  by  these  processes,  the 
core  is  heated  due  to  gravitational  contraction  and  at  a 
temperature  of  about  10   K  (corresponding  to  an  effective 

thermal  energy  of  about  200  keV) ,  helium  burning  is  inlti- 

1 P 

ated.   The  parameters  for  the  3a-reaction  leading  to   C 

are  known  to  a  fair  degree  of  accuracy,  and  hence  the 
rates  of  these  reactions  have  been  determined  (3). 

The   C(a,y)   0  reaction  in  the  helium  burning 
chain  is  believed  to  proceed  through  the  7.12  MeV  (1~) 
level  in   0  which  is  about  4-0  keV  below  the  threshold. 
Very  little  is  known  about  this  state,  especially  the 
reduced  alpha  particle  width  which  Is  necessary  in  the 

19        1  f\ 

determination  of  the   C(oc,y)   0  reaction  rate  at  stellar 


2 

temperatures.   The  reduced  alpha  particle  width  9a  in  units 

of  the  Wigner  limit  is  defined  here  by  9a2  =  Ya2(3h2/2MR  )-1. 

o 

The  quantity  ya  is  related  to  the  partial  width  Ta   by  the 

2   2  2 

expression  T     =  2y  An  ,  where  An   is  the  penetrability.   Alpha 

particle  model  calculations  (where  the  iO0  nucleus  is  assumed 
to  be  formed  by  four  alpha  particles  at  the  corners  of  a 
tetrahedron)  indicate  that  the  reduced  width  for  the  7.12  KeV 
state  is  near  unity  (*0;   Cluster  model  calculations  by  Roth 
and  Wildermuth  indicate  a  larger  variation  for  the  reduced 
alpha  width  (5).   The  value  of  the  reduced  width  for  states 
obtained  by  assuming  a  XfcC  +  He  cluster  (in  their  ground 
states)  is  about  0.6.   If  the   C  nucleus  is  in  the 
4.^3  KeV  (2+)  state  then  the  reduced  width  is  about  0.02. 
Roth  and  Wildermuth  assign  the  7.12  HeV  level  to  the  12C  +  ^He 
cluster,  indicating  that  the  reduced  width  is  nearer  the 
value  0.6.   Since  the  final  abundances  of  12C  and   0  will 
depend  critically  on  the   C(<x»y)   0  reaction  rate  (and  hence 
on  0   ) ,  a  more  accurate  determination  of  the  reduced  alpha 
particle  width  of  the  7.12  KeV  level  is  required.   It  is 
possible  to  obtain  some  information  on  9a  by  examining  the 

C(cxfy)   0  reaction  cross  section  at  alpha  particle  bom- 
barding energies  near  2  MeV,   Because  there  were  no  accurate 
estimates  of  this  reduced  alpha  width,  Delnzer  and  Salpeter 

(6),  and  also  Fowler  and  Hoyle  (7),  calculated  the  final 

12 

abundance  of  fcC  in  the  core  of  a  helium  burning  star  for 

various  assumed  values  of  9   . 
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Prior  Results 

Prior  results  on  the  12C(a,Y)l60  reaction  have  been 
limited  to  work  done  by  Allan  and  Sarma  (8),  Bloom  at  al.  (9), 
Heads  and  Mclldowie  (10),  and  Larson  and  Spear  (11,  12).   In 
the  preliminary  work  done  by  Allan  and  Sarma,  alpha  particles 
(Ea  =1.6  MeV)  were  used  to  bombard  thick  carbon  targets 
while  observing  the  high  energy  gamma  rays  with  a  5.08  cm 
thick  Nal(Tl)  crystal.   The  yields  were  determined  by  comparing 
the  counting  rates  in  the  high  energy  region  with  the  known 
counting  rate  due  to  the  h.kj,   MeV  gammas  from  the  9Be(a,nY)12C 
reaction.   They  obtained  the  value  3  x  10'5   barns  for  the 
integrated  cross  section  at  Ea  =  1.6  MeV.   It  is  now  felt  that 
the  correct  cross  section  is  much  smaller  than  this. 

Bloom  et  al.  examined  the  1"  state  In  l60  at  9.59  MeV. 
They  were  investigating  isotopic  spin  mixing  of  this  T  »  0 
state,  from  which  the  ground  state  radiation  should  be  inhib- 
ited. The  radiative  width  Ty   for  this  state  was  found  to  be 
0.006  eV  corresponding  to  |H|2(E1)  =  1  x  10-5  in  weisskopf 
units.  From  this  it  was  concluded  that  there  was  some  T  =  1 
mixing  present. 

Meads  and  Mclldowie  studied  the  2+  excited  states  in 
160  at  9.85  MeV  and  11.52  MeV  in  order  to  obtain  the  radiative 
widths  of  these  states.  Using  enriched  12c  targets  and  a 
10.2  cm  x  10.2  cm  Nal(Tl)  crystal,  the  values  0.02  +  0.01  and 
0.9  +  0.2  eV  were  obtained  for  these  radiative  widths, 
respectively.   The  state  at  9.59  MeV  was  observed,  but  the 
radiative  width  was  not  found  for  this  state. 
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Larson  and  Spear  examined  the  region  Ea  =  2.8  to  8.3  HeV 
in  some  detail  with  enriched  12C  targets  using  a  10.2  cm  by 
10.2  cm  Nal(Tl)  detector.   The  region  of  interest  for  the 
present  experiment  Includes  principally  the  resonances  at 
Ea  =  3.24  and  3.58  MeV  corresponding  to  excited  states  in 
l60  at  9.59  HeV  (1~)  and  9.85  MeV  (2+) ,  respectively.   The 
radiative  width  for  the  9.59  MeV  state  was  given  by  Larson 
as  0.022  +  0.005  eV,  and  the  value  of  the  radiative  width 
for  the  9.85  MeV  state  was  0.0059  ±  0.0006  eV.   Also,  the 
cross  section  at  Sa  =  3. 2k   HeV  was  given  as  36   nb  (10-9  barns) 
(11).   The  disagreement  of  Larson  with  the  results  of 
Bloom  et  alA  was  explained  by  assuming  the  latter  had  inad- 
vertently lost  data  due  to  their  background  subtraction  or 
renormallzatlon.   Larson  and  Spear  were  In  better  agreement 
with  the  expected  isobaric  spin  impurity  predicted  by 
Wilkinson  (13).   it  was  not  possible  for  them  to  determine 
the  capture  cross  section  for  lower  alpha  particle  bombarding 
energies.   Cascade  radiation  of  about  7  MeV  was  observed 
from  the  9.85  MeV  state;  however,  it  was  Impossible  to 
determine  whether  this  cascade  proceeded  through  the  7.12  MeV 
(1-)  state  or  the  6.92  MeV  (2+)  state  in  l60.  A  width  of 
0.0012  +  0.0004  eV  was  given  for  this  radiation,  assuming  an 
isotropic  distribution. 

A  recent  endeavor  on  the  12C(a,y)l60  reaction  was 
presented  at  the  1968  American  Physical  Society  Meeting  in 
Washington  by  Adams  et  al,,  (14).   Using  time-of -flight  methods 
to  discriminate  against  the  neutron  background,  a  value  of 
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10  nb  was  obtained  for  the  capture  cross  section  at 
Ea  =  2.75  MeV.   Data  accumulation  was  hindered  somewhat  by 
low  yields  resulting  from  the  detector  being  relatively  far 
from  the  target  and  by  low  beam  currents  from  the  tandem 
accelerator  employed  in  the  experiment. 

Intent  of  Experiment 

At  the  time  this  experiment  was  first  proposed  it  was 
hoped  that  absolute  cross  section  measurements  could  be 
extended  to  an  alpha  particle  bombarding  energy  below  1.5  MeVj 

however,  this  did  not  prove  feasible  because  of  the  extremely 

•jo       i  c. 

low  yield  from  the   C(c,y)   0  reaction,  and  the  lowest 

bombarding  energy  at  which  the  cross  section  was  determined 
was  1.9  MeV. 

.Because  of  the  discord  between  the  results  of  Larson 
and  those  of  previous  experimenters  concerning  the  radiative 
widths  of  the  9.59  MeV  and  9.85   MeV  excited  states,  it  was 
felt  that  these  radiative  widths  should  be  measured.  Also, 
since  Larson  was  not  able  to  determine  the  state  through 
which  the  cascade  radiation  from  the  9. 85  MeV  state  in  *°0 
proceeded,  it  appeared  desirable  to  examine  this  cascade 
radiation  using  a  20  cm^  Ge(Li)  detector. 


CHAPTER  I 


THEORY 


Nuclear  Processes  In  Stars 

In  search  of  an  -understanding  of  the  origin  of  the 
elements  in  the  universe,  many  theories  have  been  formulated. 
So  far  the  theories  fall  into  two  classes.   One  class  assumes 
that  the  elements  were  formed  in  a  primordial  state  of  the 
universe.   Two  examples  of  this  class  are  the  non- equilibrium 
theory  of  Alpher  and  Herman  (15)  and  the  poly-neutron  theory 
of  Mayer  and  Teller  (16).   The  other  type  proposes  that  stars 
are  the  principal  instrument  in  element  formation.  A 
relatively  complete  theory  of  this  type  has  been  formulated 
by  Burbidge,  Burbidge,  Fowler,  and  Hoyle  (17).   Primordial 
theories  distribute  the  elements  on  a  cosmic  scale;  however, 
the  theories  Imply  that  the  distribution  should  be  independent 
of  time.   This  Is  contrary  to  what  is  observed  since,  for 
example,  anomalies  (such  as  the  presence  of  technetium) 
have  been  observed  in  certain  stars.   Stars  may  distribute 
material  by  ejection  (such  as  the  explosion  of  a  supernova) 
which  would  result  in  a  distribution  of  the  elements  on  a 
cosmic  scale. 
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The  theories  must  all  be  able  to  explain  the  relative 
abundance  distribution  of  the  elements.   Suess  and  Urey  have 
constructed  a  "smoothed"  abundance  curve  by  employing  data 
obtained  from  terrestrial,  meteoritic,  and  solar  measurements 
(18).   It  is  impossible  to  say  unequivocally  that  such  a 
curve  is  universal.   Thus,  one  only  attempts  to  explain  the 
origin  of  the  material  that  is  observed.   The  general 
feature  of  the  abundance  curve  is  that  it  decreases  exponen- 
tially from  A  =  1  to  about  A  =  100  (figure  1).   This  is 
from  Burbldge  et  al.  (17).   Other  principal  features  Include 
the  change  in  slope  at  A  =  100,  the  scarcity  of  D,  Li,  Be  and 
B,  the  relatively  high  abundance  of  alpha  particle  nuclei 
such  as  °0,   Ne,   Ti ,  large  peaks  centered  at  A  equal  to 
about  85,  13^,  202,  and  finally  the  small  abundance  of 
proton-rich  heavy  nuclei.  Any  complete  theory  must  be  able 
to  explain  all  of  these  features. 

Burbidge,  Burbidge,  Fowler,  and  Hoyle  proposed  eight 
processes  in  stars  that  would  account  for  these  features  in 
the  abundance  curve  (17).  A  few  modifications  to  the  theory 
have  occured  since  the  original  article,  and  the  number  of 
processes  has  increased.   The  processes  are  labeled  hydrogen 
burning,  helium  burning,  carbon  and  oxygen  burning,  the  alpha 
process,  the  equilibrium  process,  the  s-process,  the  r-process, 
the  p-process  and  the  P-process.   It  will  be  shown  that  the 
12C(o.iY)   0  reaction  is  in  the  helium  burning  chain  of  reac- 
tions.  However,  to  better  understand  the  motivation  behind 
this  experiment  it  is  desirable  to  examine  the  over-all 
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Figure  1.   "Smoothed"  abundance  curve  of  the  elements, 
Burbidge  et  al.  (17). 
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theory  slightly.   Also,  all  of  the  separate  processes  inter- 
connect to  form  the  complete  theory.   It  is  for  these  reasons 
that  a  brief  description  of  the  individual  processes  will  be 
given. 

The  P-process  is  introduced  to  explain  the  formation 
of  the  elements  lithium,  beryllium,  and  boron.   These  elements 
are  rapidly  reduced  to  helium  through  numerous  nuclear  reac- 
tions with  protons  in  stellar  interiors;  thus,  it  appears 
that  they  must  be  produced  near  the  surfaces  of  magnetic 
stars  through  spallation  of  carbon  by  highly  energetic  protons, 

The  p-process  is  employed  to  explain  the  existence  of 
certain  proton-rich  nuclides  found  in  the  heavy  elements. 
These  proton-rich  elements  are  formed  by  (p,y)  or  [y,n) 
reactions  on  the  elements  along  the  nuclear  stability  line. 

Hydrogen  burning  is  the  synthesis  of  helium  from 
hydrogen  (table  1)  (19).   The  alternate  endings  seem  to 
depend  on  the  temperature  of  the  region  where  the  hydrogen 
burning  is  occurring.   The  quantity  SQ  is  the  cross  section 
factor  at  stellar  temperatures  (19).   The  cross  section 
factor  is  directly  proportional  to  the  reaction  rate  which 
determines  abundances.   The  cross  section  factor  S  is 
related  to  the  cross  section  a   by 

S(E)  =  o(E)Eexp(31.28  Z-^A^E*^)  keV  barns, 

where  a(E)  is  the  cross  section  in  barns  (10"2^  cm2) 
measured  at  the  center-of-mass  energy  E  in  keV  (17). 


10 


o\ 


PQ 


1 

o 

-p 
o 

u 

ft 

I 

o 
-p 
o 


Q) 
H 


I  H 


CQ 

a 

I 

i 
> 

© 


W 


CM 

CM        -3- 
I  I 

o      c 

H        r- 


H 
c^ 


M 
o 


CO 

o 

H 

H 
H 


© 


CM 
M 

on 


CM 


CN 

•3- 


VN 


ra 

CO 

o 

H 

I 

CM 


© 

S3 

^ 

Si 

© 

-P 

03 

CO 

00 

MJ 

CM 

Gj 

CO 

ON 

00 

CM 

c 

© 

• 

• 

• 

•H 

H 

CM 

o 

CM 

MD 

© 

H 

iH 

CM 

ca 

h 

a 

It 

u 

o 

>» 

I 


o 

03 

*-* 

03 

*-* 

>> 

o: 

aj 

© 

-d 

a 

o 

o 

CM 

£ 

rH 

O- 

-3"         U 

o 

cv] 

O       1  h 

i-l 

«— > 
03 
CO 

o 

H 

1 

•— * 

vtn 

^1 

O- 

O 

n 

VO 

o 

CN 

V>. 

rH 

CM 

C^N 

1 

CO 

H 

o 

H 

1 

r-A 

• 

O 

H 

H 

X 

II 

II 

O 

0^ 

|h 

It- 

O 

3 

>- 

a 

+ 

• 

* 

rH 

A 

« 

+       © 

© 

H 

■H         W 

H 

CO 

+ 

CM 

J- 

►-1   -* 

.Q 

fl 

O-          CM 

N, 

a} 

o 

+ 

>-       + 

+ 

f     t 

>- 

H 

■H 

•P 
O 

+ 

+        1> 

© 
cq 

-r 

+ 
+ 

t 

©  -3" 

c*- 

I        w 

pq 

31 

P 

f     i 

h           * 

©     H 

u       co 

^ 

CM 

O               T 

© 

0    ! 

+           © 

ft 

r 

+        + 

©         «H 

©        CM 

ft 

w 

w 

PQ      ^ 

w 

ra 

w 

W    CN. 

J- 

o-  v  o- 

rH 

CO           | 

t .' 

1        © 

© 

EH 

pH 

+ 

rH 

+  

>  + 

\ 

>    + 

© 

© 

© 

w 

Q     ^W 

M 

pq 

PQ      cq 

rH 

CM       C*^ 

cn 

o- 

CO        CO 

03 
CO 

o 

I 

On 
CM 


c«^ 

rH         O-        Ol       H 
O         O-        H       ON 


> 

© 


c^» 


MO 

CM 


(S 
■P 
O 

EH 


© 

J* 


11 

Z-j_  and  Z0  are  the  charges  of  the  Interacting  particles  In 
units  of  proton  charge,  and  A  =  A]A0(Ai  +  Aq)"1  is  the  reduced 
mass  of  the  system  in  atomic  mass  units.   S  results  after 
effects  due  to  resonances  and  barrier-penetration  have  been 
eliminated  in  the  expression  for  the  cross  section. 

In  stars  that  contain  enough  carbon,  nitrogen,  or 
oxygen,  helium  may  be  produced  from  hydrogen  through  the 
carbon-nitrogen- oxygen  cycle.   The  cycle  starts  with  radiative 
proton  capture  on  -^C  and  progresses  to  the  reaction 
^N(pta)^C.  After  the  hydrogen  in  the  core  has  been  exhausted, 
the  core  begins  to  contract  due  to  gravity.  As  this  occurs 
the  temperature  rises  from  the  10?  °K  that  was  present  during 
the  hydrogen  burning  process. 

When  the  temperature  is  about  10^  °K  the  core  begins 
to  burn  helium.   The  principal  reactions  involved  are  shown 
in  figure  2  (19).   The  0+  ground  state  In  8Be  and  the  0+ 
state  at  7.65   MeV  in  12C  allow  the  first  two  reactions  to 
proceed  at  a  reasonable  rate.   The  formation  of  *"0  appears 
to  be  caused  mainly  by  the  1"  state  at  7.12  MeV  (figure  2). 
This  is  about  *J-0  keV  below  the  12C  +  ^He  threshold.  Using 
resonance  fluorescence,  Swann  and  Metzger  found  the  mean 
lifetime  of  this  state  to  be  about  1.  0  x  10~1/f  s  (as  corrected 
from  their  previously  obtained  value),  corresponding  to  a 
radiative  width  of  about  0.065  eV  (20).   The  alpha  particle 
width  for  this  state  remained  unknown.   If  the  reaction  is 
due  to  a  single  resonance  (as  is  the  case  for  12C(<x.y)  0), 
then  it  can  be  shown  that  the  cross  section  factor  S  is 
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Figure  3.  Energy  levels  of  ^O.   Lauritsen  and 
Ajzenberg-Selove  (21). 
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proportional  to  92  (the  reduced  particle  width  of  that  state) 
(17).   Thus,  it  is  evident  that  a  determination  of  the 
reaction  rate  requires  a  determination  of  that  reduced  width. 

Figure  3  shows  the  energy  levels  of  l60  (21),   It  is 
possible  for  the  broad  (Tcla  =  6>+5   keV)  1"  level  at  9.59  MeV 
to  interfere  with  the  1-  level  at  7.12  MeV.   This  inter- 
ference will  either  enhance  or  decrease  the  capture  cross 
section  in  the  region  between  these  resonances,  depending 
on  whether  the  interference  is  constructive  or  destructive, 
respectively.   If  the  relative  sign  of  the  amplitudes  of 
these  states  is  minus,  constructive  interference  will  occur. 
If  the  sign  is  positive,  destructive  interference  will  occur. 
Thus,  a  determination  of  9a2  is  possible  by  fitting  the 
observed  cross  section  in  the  region  Ea  =  1.5  to  3  MeV  to 
calculated  values  of  the  cross  section  (for  assumed  values 
of  ea  ,  and  computing  both  destructive  and  constructive 
values),  A  calculation  of  this  type  has  been  done  by 
Tombrello  (22).   It  Is  therefore  necessary  to  measure  the 
capture  cross  section  in  this  region  in  order  to  determine 
the  reduced  alpha  particle  width  ea2.   Larson  has  measured 
the  cross  section  at  an  alpha  particle  bombarding  energy 
of  3.24-  MeV,  but  this  does  not  allow  a  determination  of 
the  interference  effects  (11). 

The  relative  abundance  of  l60  to  12C  will  depend  on 
the  cross  section  of  this  reaction  at  stellar  temperatures? 
thus,  the  necessity  for  a  better  extrapolation  forms  an 
Incentive  for  the  present  experiment.   The  main  difficulty 
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In  performing  the  experiment,  besides  the  extremely  low  cross 
section  values  for  the  ^C(aiY)  0  reaction,  is  the  relatively 
large  cross  section  (20-100  mb)  for  the  ^CCa.n)1  0  reaction, 
since  Nal(Tl)  detectors  are  quite  sensitive  to  neutrons. 

After  the  helium  has  burned,  the  core  again  contracts, 
and  the  stars  will  begin  burning  the  carbon  and  oxygen. 
In  this  way  nuclides  up  to  2°S1  may  be  formed. 

As  the  temperature  rises  (2  x  10^  °K)  photodisintegra- 
tion  occurs  producing  alpha  particles,  protons,  and  neutrons. 
Since  the  thresholds  for  (ytP)  and  (Y.n)  processes  on  nuclei 
with  A  =  2Z  =  2N  are  higher  than  for  other  nuclei,  there 
tends  to  be  an  Increase  In  abundance  for  nuclides  with 
A  =  2Z  =  2N,   These  give  the  appearance  of  being  composed 
of  alpha  particles.   This  Is  known  as  the  alpha  process. 

As  the  temperature  rises  to  about  ^  x  10^  °K,  many 
nuclear  reactions  are  occuring  so  rapidly  that  a  statis- 
tical equilibrium  is  created.   This  is  known  as  the  e-process. 
Since  electron  capture  can  occur  quickly  at  this  tempera- 
ture, there  is  a  shift  to  nuclei  that  have  N  slightly  larger 
than  Z  in  the  vicinity  of  the  "iron  peak." 

The  terms  "r-and  s-processes"  refer  to  rapid  and 
slow  neutron  capture,  respectively.   The  s-process  is  the 
capture  of  neutrons  such  that  there  Is  enough  time  between 
captures  to  permit  beta  decay.   The  necessary  neutrons  for 
the  s-process  are  produced  by  (afn)  reactions  on  nuclei  such 
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n      21 
as  JC   and   Ne.   Since  certain  elements  beyond  bismuth  have 

decay  times  that  are  short  compared  to  the  characteristic 

time  for  the  s-process,  the  elements  heavier  than  bismuth 

could  not  be  produced  by  the  s-process.   But  neutron  capture 

at  a  rapid  rate  (the  r-process)  could  produce  the  heavy  elements 

such  as  uranium,  etc.   Since  the  cross  section  in  the 

r-process  and  s-process  are  small  at  neutron  magic  numbers, 

the  abundances  are  large  at  these  points.   The  r-process 

produces  the  peaks  in  the  abundances  curve  at  lower  values  of 

A  than  the  s-process,  of  course,  since  the  r-process  occurs 

to  the  neutron  rich  side  of  the  stability  (N  =  Z)  line;  thus, 

the  formation  of  the  double  peaks  observed  in  the  abundance 

curve  may  be  explained  by  these  processes. 

As  the  Iron-group  nuclei  increase,  the  star  again 

contracts  and  the  core  implodes;  however,  nuclear  reactions 

are  still  occurring  in  the  outer  part  of  the  star.   This 

implosion  creates  the  supernova.   This  also  produces  a  rapid 

increase  in  temperature,  causing  an  emission  of  a  large 

part  of  the  star;  furthermore,  neutrons  are  produced  profusely 

and  are  thus  convenient  for  the  execution  of  the  r-process. 

In  this  way  evolution  of  the  star  is  completed. 

Radiative  Alpha  Capture 

Radiative  alpha  capture  has  often  been  employed  for 
determining  detailed  information  (such  as  the  energy,  spin, 
parity,  etc.)  on  excited  states  of  nuclei  (23,  2^-).   These 
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parameters  (along  with  others  such  as  radiative  widths, 
branching  ratios  of  gamma  rays,  mixing  parameters,  etc.) 
are  required  in  comparing  nuclear  models  and  extending 
their  usefulness.   The  alpha  particle  capture  reaction  is 
different  from  proton  capture  in  that  the  alpha  particle 
Is  spiniess;  thus,  there  is  a  simplification  in  angular 
correlation  meaurements.   If  the  target  nucleus  is  also  spln- 
less,  as  is  the  case  for  even-even  nuclei  such  as  xo0,  a 
further  simplification  occurs.   In  this  case  only  the  relative 
motion  of  the  interacting  particles  may  contribute  to  the 
spin  of  the  compound  state.   This  requires  that  the  possible 
values  of  spin  (J)  are  limited  to  J  =  1,  2,  ?>,,.. 
These  are  simply  the  eigenvalues  U)    of   the  orbital  angular 
momentum  operator.   The  concept  of  parity  is  determined  by 
the  behavior  of  the  wave  function  for  the  system  upon 
reflection  of  spatial  coordinates  (25).  If  upon  this  reflec- 
tion the  spatial  part  of  the  wave  function  does  not  undergo 
a  sign  change,  the  parity  (tt)  is  even  (+).   If  the  sign 
changes,  the  parity  is  odd  (-).   Furthermore,  It  is  known 
that  the  parity  is  determined  by  the  orbital  angular 
momentum  eigenvalue  {!)*      If  /  is  even  (0,  2,  Urt    .  .  .), 
than  the  sign  of  the  wave  function  will  not  change  upon 
spatial  reflection,  and  the  parity  will  be  even  (+).   If 
/  is  odd  (1,  3,  .  .  .),  the  sign  will  change,  and  the 
parity  will  be  odd  (-).   Hence,  in  cases  where  the  pro- 
jectile and  target  are  both  spiniess,  the  spin  (J)  and 
parity  (rr)  of  the  states  that  can  be  excited  are  limited 
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to  Jn  =  0+,  1",  2+,  3",  ,  ,  ,   These  are   the  so-called 

states  of  "natvixal"  parity. 

For  this  experiment  the  concept  of  the  compound 
nucleus  Is  also  useful  (25).   Here  the  incident  particle  is 
captured  and  a  compound  system  is  formed,  with  the  energy 
shared  among  all  the  particles  of  the  compound  system.  A 
relatively  long  time  elapses  before  sufficient  energy  is 
concentrated  onto  a  particle  to  allow  It  to  escape.   This 
time  is  long  compared  to  the  "characteristic  nuclear  time," 
which  is  the  time  required  for  the  incident  particle  to 
traverse  the  target  nucleus  without  interacting.   Since  the 
decay  time  is  long  compared  the  the  "characteristic  time," 
the  compound  nucleus  "forgets"  its  method  of  formation.   If 
there  is  sufficient  energy,  the  compound  nucleus  can  decay 
through  particle  emission  in  addition  to  gamma  emission. 
The  states  excited  in  this  way  are  called  virtual  states. 
If  there  is  not  sufficient  energy  to  allow  particle  emission, 
the  compound  nucleus  can  decay  by  gamma  emission  or  by 
internal  conversion.   In  any  case  conservation  laws  are 
obeyed. 

In  gamma  ray  emission,  the  angular  momentum  of  the 
gamma  ray  for  the  ground  state  transition  must  be  equal  to 
the  angular  momentum  of  the  compound  system,  which  is  equal 
to  the  orbital  angular  momentum  in  the  particular  case  where 
the  spin  of  the  interacting  particles  is  zero.  A  light 
quantum  has  an  angular  momentum  equal  to  hf[!   +  1)^.   The 
maximum  projection  is  fft.   The  multipole  order  of  radiation 
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is  2  ;  thus,  for  example  in  the  case  where  the  interacting 

particles  are  spinless,  p-wave  capture  will  result  In 
electric  dlpole  radiation  for  the  ground  state  transition. 
Since  the  radiation  is  transverse,  there  cannot  be  an  /  =  0 

multipole.   This  is  the  reason  that  0+ >  0+  radiative 

transitions  are  strictly  forbidden. 

Radiative  Capture  of  Alpha  Particles  by  jjp 

The  two  resonance  levels  in  160  that  can  be  reached 
by  the  present  experiment  are  the  2+  level  at  9.85  MeV  and 
the  1"  level  at  about  9.59  MeV  (figure  3).   These  occur  at 
incident  alpha  particle  energies  of  3.58  and  3.2*f  MeV, 
respectively. 

The  level  at  9.85  MeV  has  a  total  width  rcffl  of  0.75  keV 
(26,  27).   This  narrow  resonance  proved  quite  useful  for 
determining  the  target  thickness.   The  thickness  of  the 
target  is  essential  in  the  determination  of  cross  sections 
and  radiative  widths.  Also,  because  of  the  large  cross 
section  for  the  13C(a,n)l60  reaction  Immediately  preceding 
this  resonance,  the  XCC   enrichment  was  indicated  through  a 
comparison  of  the  excitation  function  for  a  natural  carbon 
target  with  that  from  an  enriched  12C  target. 

Unfortunately  a  broad  peak  due  to  the  13C(a,n)l60 
reaction  is  located  in  the  same  alpha  particle  energy  region 
as  the  9.59  MeV  peak  (28).  This  makes  the  gamma  radiation 
from  the  9.59  MeV  level  extremely  difficult  to  observe. 
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The  9.59  MeV  level  Is  a  broad  resonance  having  a  total 
width  Tcm  of  6k5   keV  (26).   Hence,  the  peak  extends  over  a 
considerable  energy  region,  but  the  determination  of  the 
off-resonance  yield  is  hampered  by  the  sensitivity  of  the 
Nal(Tl)  spectrometer  to  neutrons  originating  in  the 
-^C(a,n)l60  reaction. 

Below  about  Ea  =  2.9  MeV  the  main  process  is  nonreso- 
nant  capture.   Here  the  cross  section  for  the  12C(a,Y)l60 
reaction  is  extremely  small  (<10  nb) ,  while  the  13C(a,n)l60 
reaction  remains  relatively  large  (-20-100  tab).   Thus,  it  is 
even   more  difficult  to  determine  the  nonresonant  reaction 
cross  section  at  the  lower  energies.   Around  Ea  =  2.0  MeV 
interference  effects  between  the  7.12  MeV  and  the  9.59  MeV  1" 
states  in  1*0  become  Important.   The  cross  section  will  be 
increased  or  decreased  depending  on  whether  this  Interference 
is  constructive  or  destructive.   The  magnitude  of  this  effect 
is  dependent  upon  the  parameters  of  the  states  involved. 
Tombrello  has  calculated  this  effect  assuming  ea2  =0.1, 
however,  without  knowing  the  relative  sign  of  the  amplitudes 
it  is  impossible  to  determine  whether  this  effect  ls  con- 
structive of  destructive  (22). 

Aggnlar^lstrl  but  ions 

In  the  appendix  the  efficiency  of  the  Nal(Tl)  crystal 
is  derived  assuming  an  Isotropic  distribution  of  gamma 
radiation.   It  is  necessary  to  determine  the  effects  of  an 
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anisotropic  distribution  on  these  results.   Gamma  radiation 
resulting  from  transitions  between  natural  states,  i.e. 
0  '  1  '  2  •  •  •  ♦  wl11  have  an  angular  distribution  given 
by  W(e)  =  anPn(cose),  where  Pn(cose)  is  the  nth  order  Legendre 
polynomial,  and  e  is  the  angle  between  the  emitted  photon  and 
the  incident  particle.   Since  the  ground  state  of  l60  is  0+f 
the  observation  of  radiation  Yo  to  the  ground  state  resulting 
from  the  capture  of  an  alpha  particle  in  most  cases  Indicates 
that  the  spin  and  parity  of  the  resonance  state  is  either  1" 
or  2  .  J  =  0  would  not  decay  to  the  0+  ground  state,  and 
J  =  2  is  more  likely  to  decay  to  higher  excited  states  which 
would  then  decay  to  the  ground  state  in  160.   In  particular, 
for  0+_^  2+_^  0+  and  o+_^  l"^  0+  transitions  the 
theoretical  angular  distributions  will  be  of  the  form 

W(9)^sin2e-~l-P2(cose)  (Jnal-) 

W(e)^sln22e-l+(5/7)P2(cosG)-(12/7)P4(cose)   (jn=2+). 

Prom  figure  k   it  can  be  seen  that  these  distributions  are 
entirely  distinguishable.  The  simplicity  of  these  distri- 
butions results  from  the  fact  that  there  can  be  no  radiation 
mixing,  channel  spin  mixing,  or  orbital  angular  momentum 
mixing  for  these  transitions.   However,  the  distribution 
resulting  from  the  transition  between  the  2+  states  at 
9.85  MeV  and  6.92  MeV  in  l60  (0+  d  2+  M1E2  2+j  lg 
complicated  by  the  fact  that  Ml  -  E2  radiation  mixing  is 
Possible.   In  this  case  it  is  possible  to  obtain  from  the 
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angular  distribution  the  value  of  the  mixing  parameter  X, 

where  X  =  S(S2)/S(M1)..  S(E2)  and  S(M1)  are  the  amplitudes  of 
the  S2  and  Ml  transition  probabilities,  respectively  (29). 

In  determining  the  total  yield  (and  therefore  radiative' 
widths,  cross  sections,  and  branching  ratios)  it  is  essential 
to  have  a  knowledge  of  the  angular  distribution  of  the  gamma 
radiation.   If  a  measurement  is  taken  at  only  a  single  angle, 
then  an  accurate  determination  of  the  yield  can  be  obtained 
only  with  the  aid  of  a  prior  knowledge  of  the  angular  distri- 
bution, unless  the  detector  is  so  close  to  the  target  that 
the  distribution  appears  isotropic. 

The  observed  angular  distribution  W0-jjS(G)  is  then 
"smoothed"  by  the  large  solid  angle  subtended  by  the  detector. 
This  means  that  the  larger  the  solid  angle  subtended  (the 
larger  the  detector  or  the  nearer  the  detector  to  the  target), 
the  smaller  the  angular  variations  will  be.   In  other  words, 
the  observed  angular  distribution  will  have  a  decreased 
dependency  on  higher  order  Legendre  polynomials.   Wobs(9) 
has  the  form 

Wobs(9)  =  anbnPn(cos8),     |fcn|=l. 

The  coeffl cents  bn  are  the  attenuation  coefficients  described 
by  Rose  (30).   They  are  given  by  bn  =  Jn/JQ,  where  Jn  is  the 
nth  order  integral  defined  in  the  appendix.   These  Integrals 
were  evaluated  numerically  through  the  use  of  a  subprogram 
(forming  a  part  of  a  program  that  determines  angular  cor- 
relation parameters)  (31)  (table  2).   In  general,  the  values 
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of  Jn  and  bn  decrease  rapidly  with  increasing  n  for  large 

solid  angles.   It  is  obvious  that  the  proximity  of  the  detector 
to  the  target  determines  how  rapidly  the  bn  approach  zero 
with  increasing  n.   If  the  detector  is  very  close  to  the 
target,  the  observed  distribution  will  be  quite  independent 
of  higher  order  terms. 

Absolute  Yields 

The  absolute  total  yield  of  a  nuclear  reaction  is 
defined  here  as  the  number  of  disintegrations  per  incident 
particle.   In  the  following  discussion  the  notation  adheres 
to  that  used  in  the  appendix.   The  yield  may  be  related  to 
Nobs^9^  (the  number  of  counts  in  the  spectrum  peak  observed 
at  the  laboratory  angle  e  per  incident  particle,  corrected 
for  background  and  dead-time)  by  employing  the  ideas  derived 
above  and  in  the  appendix.   The  number  of  counts  in  the 
peak  (here  defined  as  the  region  (Ex  -  1.02  MeV)  to  1.1BX) 
at  angle  e  per  steradian  per  incident  particle,  is  given  by 

N0bs(9)/^r  =  Hp  x  (YAti)  w0bs(e)/ao, 

where  r?p  =  B»?t  =  R(f>r)  -  R(J0/2) 

R   =  peak-to-total  ratio 

f   =  OAtt 
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£1  =  solid  angle  subtended  by  detector 

^  a  source  full  energy  peak  efficiency 

>7t  -  source  intrinsic  efficiency 

»7T  =  source  incident  intrinsic  efficiency 

JQ  =  the  zeroth-order  Integral  as  defined  in  the 
appendix. 

The  absolute  total  yield  is  denoted  by  Y;  hence,  Y/4tt  is  the 
average  yield  per  steradian.  Thus,  from  the  above  equation, 
it  is  seen  that  the  yield  is 

Y  =  Nobs(e)(H^Jo/2)  ^Kan/a0)(Jn/Jo)Pn(cose))-1, 

The  expression  Nobs(e)/4n  may  be  defined  as  the  "ob- 
served differential"  yield  and  denoted  by  nobs(e).   Similarly, 
the  expression  U(e)An  may  be  defined  as  the  actual  differ- 
ential yield  and  set  equal  to  n(e).   Using  the  expression 
derived  for  the  yield  it  is  easy  to  show  that  the  relation 
between  n(e)  and  nQbs(e)  is 


n(e)  = 


Gobs(e)    S(an/a0)Pn 


R(Jo/2)    S(an/a0)(Jn/J0)Pn  ' 


2? 

Radiative  V/iclth  Determination 

For  a  resonant  radiative  capture  reaction  the  cross 
section  a  is  assumed  to  depend  on  the  energy  according  to 
the  single-level  dispersion  formula  (25).   Explicitly  for  an 
(ciy)  reaction  this  has  the  form 

~  n*2cjrY  r« 


(E  -  ER)2+  r2A 
where  o>  =  statistical  weight  =  (2J  +  1)/((2I  +  i)(2S  +  1) ) 

J  =  total  angular  momentum  of  system  (spin  of  state) 

I  =  spin  of  target  nucleus 

3  =  spin  of  incident  particle 

r  =  total  width  of  resonance 


TY  =  radiative  width 


r  =  particle  width 


X  =  h(2mE)"2  =  wavelenght  of  interacting  particles 


m  =  reduced  mass  of  system  =  UaAT/(Aa  +  AT)  )1.66xl<r2/fg 
Aa  =  Incident  particle  mass  (atomic  mass  units) 
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Atp  =  target  mass  (atomic  mass  units) 

E  =  center-of-mass  energy  =  C^r/(Aa  +  A?))  x  EiaD 

Elab  -  incident  particle  energy  in  laboratory. 

At  resonance,  a  is  given  by 


aR  =  tt  x2a>ra  rY/(  r2A). 


Hence,  a  may  be  expressed  in  the  form 


aB(I,2A) 


a  = 


(e  -  eh)2  +  r2A  " 


The  relation  between  the  yield  and  the  cross  section,  assuming 
the  Breit-Wigner  formula  is  applicable,  has  often  been  cited 
in  the  literature  (32,  33,  Jk,    35).   The  relation  is  given  by 

fE 
Y  =     (a/£)dE, 

J  E-6 

where  6  is  the  target  thickness  in  energy  units  and 

£  is  the  stopping  cross  section  for  the  incident 
particle  per  dislntegrable  nucleus  in  the  target 
material. 
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The  relationship  between  £  and  the  stopping  power  Is 

£  =  l/n(dE/dx)  , 

where  -(dE/dx)  =  the  stopping  power  of  the  target  material 

-1  ^ 

n  =  gpN.M   =  number  of  disintegrable  nuclei  per  cm-7  of 

target  nuclei 

p  =  the  density  of  target  material  in  grams  per  cm^ 

N^  =  Avogadro's  number 

M  =  the  molecular  mass  of  target  material 

g  =  the  fraction  of  disintegrable  nuclei  in  target 
material. 

If  it  is  assumed  that  6  and  £  are  independent  of  E  over  the 
resonance,  then  6  =  (dE/dx)t  =  n£t,  where  t  is  the  target 
thickness  in  centimeters.   The  integral  for  Y  may  now  be 
evaluated  analytically.   The  result  is 


2TT*2o>r  r 

Y  =  SLY  x 


(E-Eg)  (E- 

arctan  -  arctan 

r/2 


E-E^n 

r/2    J' 


Tot? 

or         Y  r-     r. 


2<f 
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(E-Eg)  (E-ER-6) 

arctan  -  arctan 


r/2  r/2 


It  can  be  shown  that  this  expression  has  an  observed  width 
T '  given  by 

r  =     (r2  +  62)*. 

Thus,  the  thickness  of  the  target  tends  to  broaden  the 
observed  width.   When  the  target  thickness  is  determined 
(as  is  often  the  case)  by  examining  the  excitation  function 
across  a  narrow  resonance,  errors  enter  through  such  factors 
as  accelerator  instability,  energy-defining  slits,  etc.  The 
final  observed  width  will  thus  be  the  square  root  of  the 
sums  of  the  squares  of  all  of  these  contributions.   It  is 
evident  from  the  above  expression  for  the  yield  that  a 
maximum  occurs  at  the  energy  E  =  ER  +  6/2;  thus,  the  observed 
resonance  energy  is  shifted  to  the  higher  energy  side  by 
approximately  6/2.   This  is  only  approximate  because  in  the 
derivation  it  was  assumed  that  the  parameters  in  the  formula 
for  the  cross  section  were  independent  of  energy  over  the 
energy  integration.   This  shift  in  resonance  energy  is  one 
of  the  reasons  why  it  is  very  desirable  to  have  an  extremely 
thin  target  for  resonance  spectroscopy  measurements;  however, 
there  are  situations  where  it  may  be  more  desirable  to  have 
a  thick  target  (for  example,  to  increase  data  aquisition  rate). 
The  maximum  yield  is  given  by 

Ymax  =  (or  r/^)arctan(6/  r ). 
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In  terms  of  measured  quantities  the  resonance  cross  section 
ag  has  the  form 

£        %aX(e) 

aR  = x  - 


arctan(6/r)   R(J0/2)   (an/aQ) (Jn/J0)Pn(cose)' 

where  N^^e)  is  the  maximum  counts  observed  in  the  peak 
at  the  laboratory  angle  e.   Using  the  expression  relating 
aR  with  the  radiative  width  Ty,   it  can  be  shown  that  TY  is 
given  by 


£r         Nmax(e) 
Y 


n,  = x 


^TT^6>raarctan(6/r)   R(J0/2)  (an/aD)  (Jn/J0)Pn(cose) 

This  is  the  expression  that  was  used  to  obtain  the  radiative 
widths  from  the  observed  counting  rates.   For  off -resonance 
cross  section  measurements,  a  and  £  are  assumed  constant  in 
the  integral  over  the  energy  interval  (E  -  5)  to  E. 

Nuclear  Reaction  Rates  in  Stars 

At  a  temperature  T  the  number  of  collisions  (coll  cnT^s"1) 
between  nuclei  of  types  0  and  1,  having  a  center- of -mass 
energy  between  E  and  E  +  dE  is  proportional  to 
n0niE3exp(-E/kT)dEf  where  nQ  and  ^   are  the  number  densities 
of  particles  0  and  1,  respectively.   The  mean  reaction  rate 
(reactions  per  cubic  centimeter  per  second)  of  a  thermonuclear 
process  is  thus  proportional  to  the  integral  (36). 
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fn   n-,E2exp(-E/kT)P(E)S(E)dE, 

where  E  Is  the  center-of-mass  energy,  E2exp(-E/kT)  is  the 
Boltzmann  factor  (probability  for  the  interaction  energy  E) , 
P(E)  is  the  barrier  penetration  factor,  and  S(E)  is  the 
cross  section  factor  (defined  previously  in  section  1).   The 
quantity  S  is  a  slowly  varying  function  of  the  energy  E  for 
off -resonant  reactions.   It  can  be  shown  that  the  product 
of  exp(-E/kT)  and  P(E)  has  a  fairly  strong  maximum  centered 
at  an  energy  E0  given  by. 

E0  =  26.29(Z12ZQ2AT82)  3  keV, 

where  Z-^  =  charge  (in  units  of  proton  charge)  of  interacting 
particle  1 

ZQ  =  charge  of  particle  0 

A  =  AiAc/(Al  +  Ao)  =  reduced  atomic  mass  number  of 
interacting  particles 

Tg=  stellar  temperature  in  units  of  108  °K. 

The  quantity  EQ  Is  the  effective  thermal  energy  used  by 
Burbidge  et  al.  (17).  For  reactions  such  as  those  normally 
found  in  the  helium  burning  process,  the  effective  thermal 
energy  is  much  larger  than  the  mean  kinetic  energy. 
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This  means  that  the  region  of  interest  occurs  in  the  high 
energy  tail  of  the  Maxwell-Boltzmann  distribution  of  velocities. 

The  number  densities  n0  and  ni  are  related  to  the 
abundances  XQ  and  Xj,  expressed  as  the  fractional  amounts 
(by  weight)  of  the  interacting  nuclei t 

ni  -  NAxi  p  Ai  » 

where  FA  is  Avogadro's  number,  p  Is  the  density  in  g  cm~3,  and 
Ai  is  the  atomic  mass  number  of  species  1.   The  mean  reaction 
rate  Pj_(0)  per  nucleus  of  type  0  (expressed  in  terms  of  SQ) 
for  the  interaction  with  nuclei  of  type  1  is  explicitely 
given  by 


p1(0)  =  1/t1(0)  =  3.63  x  107SQf 


7< 


AT^J  "a^ 


e 


where  3  =  9.15(Z12Z()2A/T8)3  . 


The  quantity  t1(o)  in  the  above  expression  is  the  mean 
lifetime  (seconds)  cf  the  nuclei  of  type  0  for  interaction 
with  nuclei  of  species  1.   The  quantity  f  is  the  correction 
due  to  electron  screening  of  the  bare  nuclei  as  explained 
by  Salpeter  (37).   This  is  a  result  of  the  spherical 
polarization  of  the  electron  gas  immediately  surrounding 
the  bare  nucleus. 

If  the  reaction  rate  is  due  principally  to  a  single 
resonance  at  Er  (not  in  the  range  of  the  effective  thermal 
energy),  then  an  expression  for  S0  can  be  obtained  by  using 


3^ 
the  Breit-Wigner  single  level  formula  for  a.      The  result 
is  (17) 

2 
3Q  =  3.10  x  10-'  9    °  2  H 


A  [K2/+  ^z)]2     (Eo-Er)2+r2A 


where  9i   =  reduced  particle  width 
z  =  2(BC/EH)2 


2 


Bc  =  Coulomb  barrier  height  =  Z]Z0e  /R 
ER  =  ft2/2MH2 

K2/  +  1^  ls  the  modif,le<i  Bessel  function  of  order 
2/  +  1. 

The  quantity  a/,  is  the  function  of  z  defined  by  Burbidge 

et  al.  (17).   Values  for  the  modified  Bessel  function  have 
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been  tabulated  (38).   For  the   C(a,Y)   °  reaction  (assuming 

the  1"  state  at  7.12  MeV  in   0)  the  expression  for  SQ  is 

1.17  x  105  ea  (7.12) 
So  =  keV  barns, 

(5T8T  +  l)2 

where  the  following  values  have  been  usedi 


r2A  «  (e0  -  Er)2 


r2  =  rY   =  0.065  eV 
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Er  =  ^0  keV 


E0  =  200  T8  3" 


The  reaction  rate  is  given  by 


Pa(12c)  =  (ta(12c))-1 


1.69  x  l012xaea2f12  p        9ni.X,   »i 
=  ; — bi exp(-69.2  T8  3  )   s  •L. 


T8  3  (5T83  +  i)2 


CHAPTER  II 

EXPERIMENTAL  METHOD 

Energy  Analyzing  System 

Singly  charged  helium  ions  at  energies  up  to  k   MeV 
were  available  from  the  University  of  Florida  Van  de  Graaff 
accelerator.   Beam  currents  for  the  present  experiment  were 
between  1.5  and  3  /*A.   Some  preliminary  experiments  were 
carried  out  using  beam  currents  in  the  range  of  4-7  /xA; 
however,  the  rate  of  target  deterioration  was  Intolerable 
with  these  currents.   The  singly  ionized  He  beam  was 
deflected  through  a  90°  arc  by  the  use  of  a  momentum  analyz- 
ing magnet  (of  radius  45.72  cm).   The  magnetic  field  B 
determined  the  energy  of  the  He  beam  traversing  the  slit 
system  of  the  analyzer.   This  magnetic  field  was'  determined 
by  the  usual  nuclear  magnetic  resonance  (NMR)  methods  (23i  39) 
The  NMR  fluxmeter  (Varian  Model  F-8A)  was  supplied  by 
Varian  Associates  of  Palo  Alto,  California.   Its  upper  limit 
of  measurement  is  52  kG.   This  easily  surpasses  the  upper 
limit  requirements  of  12.6  kG  for  4  MeV  alpha  particles. 

The  NMR  unit  was  calibrated  using  the  2.4374  MeV  and 

P/i       pfi 
3.1998  MeV  resonances  in  the   Mg(a,y)   Si  reaction  (40,  4l). 
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7       7 

The  calibration  has  also  been  examined  using  the  Ll(p,n)  Be 

threshold  at  1.8807  lieV.      The  calibrations  were  accurate  to 
within  1  keV  ('!2).  A  Fortran  program  for  tabulating  NMH 
frequencies  versus  alpha  particle  energies  in  5  keV  intervals 
was  employed  in  this  experiment  (39).  Energy  variations 
due  to  the  finite  apertures  of  slits  introduces  some  error 
into  the  determination  of  alpha  particle  energies.   It  was 
assumed  that  this  error  was  +  0.2  per  cent.   The  magnitude 
of  the  error  is  not  important  for  off -resonance  yields 
(since  thick  targets  were  used)  and  enters  mainly  through 
the  determination  of  the  target  thickness  using  the  9.85   MeV 
level  in  ~60.   Current  integration  was  accomplished  through 
the  use  of  an  Elcor  current  integrator  (Model  A30QB)  (calibrated 
for  an  accuracy  of  1  per  cent  or  better).  The  target  was 
cooled  with  "refrigerated"  water  that  was  demineralized, 
thereby  minimizing  the  possibility  of  current  leakage  to 
ground.   This  precaution  appeared  desirable  due  to  the  high 
mineral  content  of  the  water  normally  supplied  to  the 
laboratory. 

Cold  Trap  and  Target  Assembly 

Initial  experiments  with  natural  carbon  targets 
were  performed  in  order  to  test  target  preparation  tech- 
niques. Data  were  also  later  compared  with  enriched   C 
spectra.   The  initial  target  chamber  was  a  9.1  cm  Inside 
diameter  (I.  D.)  glass  "T,"  separated  from  the  main  vacuum 


38 
system  by  an  in-line  liquid  nitrogen  cold  trap.   This 
apparatus  has  already  been  described  in  some  detail  by 
Bruton  (39).   It  essentially  consisted  of  a  holder  (cooled 
by  water)  with  accommodations  for  3  targets  plus  a  quartz 
viewer.   It  was  possible  to  slide  the  holder  (on  Vlton  0-rlngs) 
in  a  vertical  direction  (permitting  the  positioning  of  targets 
without  opening  the  system  to  atmospheric  pressure).   Pressures 
were  around  10"  torr  on  the  oil  diffusion  pump  side  of  the 
in-line  cold  trap.   "Teflon"  gaskets  were  used  throughout 
except  on  the  target  holder  and  in-line  cold  traps,  where 
Viton  0-rings  were  employed.   Carbon  deposition  from  organic 
vapors  was  not  apparent  with  this  system,  but  several  improve- 
ments were  devised  to  insure  better  results  with  the  enriched 
12C  targets. 

The  improved  in-line  cold  trap  (designed  by  Dr.  F.  E. 
Dunnam  and  Dr.  H.  A.  Van  Rlnsvelt)  had  several  features 
which  would  minimize  contamination  of  the  target  region  by 
organic  vapors  (figure  5).   This  consisted  of  an  in-line 
liquid  nitrogen  cold  trap  from  Sulfrian  Cryogenics,  Inc. 
(Model  338-1) ,  to  which  a  collimating  system  and  a  J6   cm 
long  polished  copper  tube  have  been  added.   The  collimator 
holders  were  mounted  on  the  inside  of  copper  couplers  that 
connected  the  cold  trap  to  the  vacuum  system.  Except  for 
the  Viton  0-rings  on  the  in-line  cold  trap  the  vacuum 
gaskets  for  the  beam  transport  system  were  of  indium,  which 
is  inherently  free  of  organic  contamination.   The  long 
copper  tube  (2,8  cm  outside  diameter)  was  in  good  thermal 
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contact  with  the  flanges  on  the  cold  trap.   The  colligating 

system  consisted  of  thr^e  20  nil  Ta  discs,  one  on  the  diffusion 
pump  side  of  the  cold  trap  and  two  on  the  target  side.   The 
initial  colligating  aperture  was  0,5  cm  in  diameter,  followed 
by  0.^  cm  and  0.3  cm  apertures  on  the  target  end.   This  system 
made  it  impossible  for  the  beam  to  strike  the  long  copper 
tube;  furthermore,  the  maximum  diameter  beam  spot  on  the 
target  was  limited  to  about  0.3  cm  by  the  collimators. 

The  target  holder  was  designed  by  Dr.  H.  A.  Van  Rinsvelt 
and  is  based  on  one  used  at  Utrecht  (figure  6),   "Refrigerated" 
itfater  flowing  between  the  outer  and  inner  copper  cylinders 
cooled  the  copper  endplate  (0.02  cm  thick)  which  then  cooled 
the  target  backing  by  conduction.   The  endplate  was  at  $5° 
with  the  beam  axis.   The  outer  cylinder  was  ^.8  cm  In 
diameter,  the  inner  was  2.15  cm.  A  Ta  Insert  fitted  into 
the  inner  copper  sleeve.  This  insert  had  a  circular  hole 
in  one  end  that  was  slightly  smaller  in  diameter  (1.77  cm) 
than  the  target  backing.  Extremely  thin  Ta  fittings  on  this 
end  held  the  target  in  position.   This  insert  (and  therefore 
the  target  backing)  was  pressed  against  the  copper  endplate 
by  tightening  the  bushing  (lined  with  Ta  on  the  beam  side). 
A  greaseless  Viton  0-ring  sealed  the  target  holder  to  the 
glass  reducer.  An  annulus  of  thin  "Teflon"  was  placed 
between  the  glass  and  the  inner  surface  of  the  flange.   This 
limited  the  flow  of  any  organic  vapors  (originating  from  the 
Viton  0-ring)  into  the  target  region.   The  cold  trap  and 
target  holder  assembly  have  the  following  characteristics t 
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a)  The  beam  traverses  a  relatively  long  path  through 
an  effective  cold  trap  "before  reaching  the  target 

chamber, 

b)  The  target  holder  is  easily  and  quickly  removed 
from  the  system. 

c)  The  target  backings  are  effectively  cooled  to 
Insure  maximum  target  life  with  relatively 
large  beam  currents. 

d)  The  detector  may  be  located  close  to  the  target 
in  taking  excitation  runs  (about  0,1  era), 

e)  It  is  possible  to  outgas  the  interior  of  the  target 
chamber  by  heating  the  surfaces  prior  to  the 
performing  of  an  experiment. 

This  assembly  effectively  minimized  carbon  deposition  on  the 

target,   Fressures  on  the  diffusion  pump  side  were  typically 

2-4  x  10   torr.   The  enriched  *2C  data  were  acquired  with 

this  apparatus. 


Detectors  and  Circuitry 

The  Nal(Tl)  crystal  (12.7  cm  x  12.7  cm)  was  enclosed 
in  a  28  cm  diameter  lead  shield  in  order  to  minimize  room 
background  (figure  7).   This  was  mounted  on  a  table  that 
allowed  the  detector  to  be  rotated  about  the  target  in  a 
horizontal  plane.   When  the  Ge(Li)  and  Nal(Tl)  detectors 
were  used  simultaneously  this  rotation  was  not  allowed.   In 
this  case  the  detectors  were  set  at  fixed  angles  of  55° 
with  the  beam  axis. 

The  detector  used  for  most  of  the  measurements  was 
a  12.7  cm  x  12.7  cm  Nal(Tl)  crystal  mounted  on  an  RCA-8055 
phot omul tipller  tube.   A  20,1  crn^  Ge(Li)  detector  (Ortec 
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Figure  7.  Nal(Tl)  crystal  and  lead  shield, 


Model  6101-20)  was  used  in  determining  the  7  MeV  cascade 
radiation.   The  resolution  of  the  Kal(Tl)  spectrometer  was 
12  per  cent  for  the  1,33  Me?  gamma  rays  from  °°Co.  For  this 
experiment,  resolution  is  defined  as  the  full  width  (in  KeV) 
at  half -maximum  divided  by  the  energy  (in  HeV)  of  the  peak 
position.   The  Ge(Li)  detector  had  a  resolution  of  0.5  per 
cent  for  this  gamma  ray. 

The  Nal(Tl)  spectra  were  measured  with  a  Nuclear 
Data  512  channel  pulse  height  analyzer  (Series  130)  (figure  8). 
The  signal  from  the  photomultiplier  went  through  a  pre- 
amplifier and  then  to  an  Ortec  Selectable  Active  Filter 
Amplifier  (Model  WtO) .     The  first  two  single  channel  analyzers 
(Ortec  Model  4-13)  were  set  to  sum  over  the  9  MeV  and  7  MeV 
radiation.  The  third  single  channel  analyzer  was  used 
(in* conjunction  with  the  gate  generator  (Ortec  Model  bo6) 
and  the  linear  gate  (Ortec  Model  *K)9))  to  eliminate  the 
lower  portion  of  the  spectrum.  In  this  way  dead  time 
corrections  (from  the  multi-channel  analyzer)  to  the  spectra 
were  less  than  two  per  cent.  A  pulse  stretcher  (inserted 
after  the  linear  gate)  was  employed  on  some  of  the  runs. 
The  calibration  and  linearity  of  the  analyzer  was  checked 
often  during  the  experiment  using  a  pulser  and  several 
known  gamma  rays.   The  following  gamma  rays  were  employed 
for  this  purpose « 

a)  0.662  MeV  photopeak  from  -^^cg 

b)  1.17  and  1.13  MeV  photopeaks  from  6°Co 

c)  2.50  MeV  sum  peak  from  6°C0 
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d)  2.75  MeV  photopeak  from  2\>a 

e)  2.225  and  hAj>Z   MeV  peaks  from  a  Pu-Be  neutron 
source  and  its  paraffin  shield 

f)  8.66  MeV  (9.17  KeV-e)  escape  oeak  from  the 
13      14  *   " 

C(p,y)  IT  reaction 

g)  9.34  MeV  (9.85  MeV-e)  escape  peak  from  the 
12      16 

C(a.iY)   °  reaction. 

The  base  of  the  phot omul tipller  tube  was  modified  to 
obtain  a  fast  signal  (from  a  circuit  obtained  from  Dr.  H,  A. 
Blue).  A  negative  fast  signal  (rise  time  of  about  40  ns) 
was  obtained  from  the  anode  of  the  8055  tube.   The  slow 
signal  was  taken  from  dynode  9  and  therefore  was  positive. 
The  purpose  of  the  fast  signal  was  to  eliminate  electronic 
pile-up  problems  using  an  EG  &  G  Pile-Up  Gate  (Model  GP100). 
The  fast  circuit  components  (all  EG  &  G)  were  set-up  to 
eliminate  all  pulses  coming  within  2  (is  of  one  another 
(figure  9).  Although  this  circuit  was  employed  for  most 
of  the  measurements,  the  count-rate  was  small,  and  thus  this 
circuit  was  not  essential  for  the  execution  of  this  experiment. 

The  signal  from  the  Ortec  Ge(Li)  detector  (20.1  cm3 
active  volume)  went  into  its  matched  Ortec  Model  118A  preamp- 
lifier.  The  signal  was  then  amplified  by  an  Ortec  Selectable 
Active  Filter  Amplifier  (Model  440)  and  processed  by  1024 
channels  of  a  Technical  Measurement  Corporation  4096  channel 
analyzer  (figure  10).   In  addition  to  the  calibration  points 
used  for  the  Nal(Tl)  spectrometer  the  Ge(Li)  spectrometer  was 
also  calibrated  using  the  second  escape  peaks  of  gamma  rays 
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^9 
from  the  reaction   Mg(a.y)  SL  The  following  gamma  rays 
were  observed: 

11.72  MeV  (12.7^  KeV(Yo)-2e) 
9.95  MeV  (10.97  MeV(Yl)-2e). 
The  ADC  of  the  102^  channel  analyzer  was  not  linear  for 
the  larger  pulses;  however,  it  did  appear  to  be  linear  up 
to  at  least  10  MeV,  Another  problem  was  a  slight  drift  with 
time  in  the  gain  of  the  ADC  of  the  analyzer.   This  resulted 
in  a  slight  broadening  of  peaks  in  the  spectra  obtained  from 
long  runs.   The  results  should  still  be  valid,  even  though 
these  problems  were  present. 

Background  Radiation 

Although  there  are  obvious  advantages  in  using  the 
spinless  alpha  particle  as  the  projectile,  there  still 
exists  an  important  disadvantage.   This  disadvantage  is 
the  prolific  neutron  background  produced  by  the  bombardment 
of  1>c  by  alpha  particles.   The  carbon  is  usually  in  the 
form  of  organic  deposits  on  the  vacuum  system,  walls,  on 
slits,  on  collimators,  and  on  the  target  itself.  Neutrons 
produced  by  the  ^Cfa.n)   0  reaction  are  captured  by  the 
Nal(Tl)  crystal,  resulting  in  an  exponentially  decreasing 
background  spectrum  extending  to  about  Ey  =  9  MeV.   The 
neutron  background  from  the  Nal(Tl)  detector  results  from 
the  radiative  capture  (by  such  reactions  as  2^Na(n,y)  Na 
and  127Kn,Y)128I)  in  the  crystal  (1).   This  is  probably 
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the  principal  difficulty  in  performing  an  (cx,y)  experiment 
using  a  Nal(Tl)  detector.  The  cross  section  for  the  -^C  +  a 
reaction  in  the  region  Ea  =  2  to  5  MeV  ranges  between  20  and 
100  mb  (4-3,  W,  h5,   46,  ^7,   ^8,  28).  This  is  indeed  an 
extremely  large  cross  section  compared  to  the  12C  +  ^Ee 
capture  cross  section,  and  it  was  thus  necessary  to  avoid 

^C  contamination.   Methods  employed  Included  the  preparation 
of  enriched  (-99.^  per  cent)  12C  targets  and  minimizing 
natural  carbon  contamination  of  the  system. 

Room  background  is  also  a  problem.  For  resonance 
measurements  the  room  background  was  not  too  Important » 
however,  when  off -resonant  values  were  obtained,  it  increased 
in  significance.   Shielding  the  detector  or  using  an  anti- 
coincidence annulus  will  minimize  this  source  of  background 
(^9*).  A  plastic  anti-coincidence  annulus  surrounding  a 
12.7  cm  x  12.7  cm  Nal(Tl)  detector  had  been  constructed  at 
this  laboratory  by  D.  E.  Wulfinghoff ,  but  was  not  operational 
when  this  experiment  was  performed.   Lead  shielding  was 
provided  for  the  detector,  and  it  was  felt  that  room  back- 
ground did  not  negate  the  results  of  this  experiment. 

Target  Preparation 

Gamma  rays  experience  little  attenuation  in  passing 
through  matter;  therefore,  self -supported  targets  are  unneces- 
sary in  gamma  spectroscopy.   This  is  a  considerable  advantage 
in  target  preparation  and  handling.   However,  because  of  the 
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prolific  neutron  background  resulting  from  the  13c(a,n)1  0 
reaction  it  was  necessary  to  make  targets  enriched  in 
12C.   The  cost  of  enriched  12C  isotope  limits  the  total  amount 
available;  thus,  methods  of  preparation  were  desired  which 
would  require  a  very  small  quantity  (of  the  order  of  10-20  mg) 
of  enriched  12C  isotope.   Two  methods  were  investigated. 
These  were  the  thermal  decomposition  and  gaseous  (high- 
frequency)  discharge  methods.   In  the  thermal  decomposition 
process  (50,   51)  the  target  backing  (10  mil  Ta)  is  heated  by 
a  flow  of  current  through  it.   When  the  foil  has  heated  and 
the  pressure  reaches  a  steady  value  of  about  5  x  10"  torr, 
the  system  is  Isolated  and  methyl  iodide  (CHoI),  or  acety- 
lene (C2H2)f  is  allowed  to  enter  the  system  until  the  pres- 
sure is  about  50  torr.   The  carbon  forms  as  a  thin,  uniform 
film  on  the  Ta  backing.   The  thickness  of  the  film  is  a 
function  of  the  decomposition  time  and  CE^I  (C2H2)  pressure. 
The  temperature  of  decomposition  is  quite  different  for  the 
two  gases.   The  methyl  iodide  cracks  readily  when  the  foil 
is  heated  to  a  very  dull  red,  while  the  acetylene  requires 
a  bright  yellow  glow  of  the  foil  to  crack  properly.   The 
latter  high  temperature  causes  the  Ta  backing  to  become 
brittle,  and  the  carbon  film  tends  to  peel  from  the  backing 
unless  handled  carefully.   Targets  thus  could  be  prepared 
using  acetylene,  but  extreme  caution  had  to  be  observed. 
Targets  could  be  easily  prepared  with  the  methyl  iodide  by 
this  method,  and  many  of  these  (along  with  a  few  prepared 
from  thermally  cracked  acetylene)  were  used  in  preliminary 
experiments  on  natural  carbon. 
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It  later  became  apparent  that  acetylene  would  be 
easier  to  synthesize  than  methyl  iodide  (52).   The  gaseous 
discharge  method  allows  the  efficient  use  of  acetylene  in 
carbon  target  preparation  (53).   This  method  was  finally 
chosen  in  preparing  the  enriched  12C  targets  and,  for  this 
reason,  will  be  discussed  in  detail.   Many  of  the  techniques 
involved  (such  as  the  preparation  of  carbides,  the  synthesis 
of  acetylene  and  the  high-frequency  discharge  method)  were 
suggested  by  Dr.  J.  L.  Duggan  of  Oak  Ridge  National  Labora- 
tory during  a  visit  (supported  by  Oak  Ridge  Associated 
Universities,  Inc.)  to  the  laboratory  by  the  author. 
The  target  preparation  method  is  conveniently 
subdivided  into  four  separate  steps  which  may  be  labeled* 

1)  synthesis  of  barium  carbonate  (BaCOo) 

2)  synthesis  of  the  carbide 

3)  synthesis  of  acetylene 

40  and  the  cracking  of  acetylene  onto  the  target 
backing.  ° 

Synthesis  of  Barium  Carbonate 

The  enriched  12C  (99. <*  per  cent)  used  in  the  prepara- 
tion of  the  targets  was  obtained  from  Oak  Ridge  National 
Laboratory  in  the  form  of  graphite.   The  apparatus  involved 
in  the  production  of  BaC03  from  elemental  carbon  is  shown 
in  figure  11.   This  mainly  consists  of  a  quartz  combustion 
tube  coupled  to  a  fretted  glass  bubbling  tube.   The  carbon 
(about  20  mg)  is  contained  in  an  Al-foil  capsule  placed  in 
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Figure  11.   Barium  carbonate  apparatus. 
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the  "plugged  Teflon"  stopcock  at  the  top  of  the  system. 
Ta  wire  mesh  in  the  center  of  the  combustion  tube  holds 
the  Al  capsule  after  it  is  released  from  the  stopcock. 
Oxygen  flows  into  the  system  through  the  glass  "T."  The 
carbon  is  heated  forming  C02  that  then  reacts  readily  with 
the  Ba(0H)2  contained  in  the  test  tube  at  the  bottom  of 
the  apparatus.   The  barium  carbonate  is  formed  as  a  white 
precipitate  in  this  test  tube.   This  is  then  recovered  by 
vacuum  filtration  and  dried.   The  efficiency  (by  weight) 
of  this  process  was  typically  greater  than  70  per  cent. 

Synthesis  of  the  Carbide 

The  carbide  is  prepared  by  placing  equal  amounts  by 
volume  of  the  barium  carbonate  and  clean  calcium  filings  In 
a  small  quartz  test  tube.   Glass  wool  is  positioned  above 
the  mixture  in  order  to  confine  the  contents  during  the 
reaction.   The  quartz  tube  (containing  the  BaCO^-Ca  mixture) 
is  connected  to  a  vacuum  pump  and  heated  slightly  with  a 
Bunsen  burner.   The  system  is  then  isolated  and  heated 
vigorously  with  an  oxygen-acetylene  flame.   The  carbide 
forms  a  black  deposit  inside  the  quartz  test  tube. 

Synthesis  of  Acetylene 

Acetylene  is  produced  by  the  addition  of  water  to 
the  carbide  (5*0.   This  process  is  carried  out  employing 
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the  apparatus  shown  in  figure  12.   This  consisted  of  a 
water  reservoir  (upper  beaker),  a  reaction  chamber  and  a  series 
of  cold  traps.   The  first  two  cold  traps  (A  and  B)  were 
alcohol-dry  ice  traps.   The  purpose  of  these  traps  13 
to  Halt  the  flow  of  water  vapor  from  the  reaction  chamber. 
Acetylene  (subliming  at  -83°C)  would  not  be  removed  by  these 
traps  (having  a  temperature  of  -78  C);  however,  it  would  be 
quickly  recovered  by  trap  C  (a  liquid  nitrogen  trap  at  about 
-196°C)  {S5)>      The  final  liquid  nitrogen  trap  D  served  to 
limit  the  possibility  of  contamination  from  the  vacuum  pump, 

The  procedure  consists  of   breaking  the  quartz  tube 
containing  the  carbide  (obtained  from  the  previous  process) 
and  immediately  placing  the  pieces  in  the  reaction  chamber 
which  is  then  evacuated  and  isolated  from  the  vacuum  pump 
by  closing  the  stopcock  farthest  from  the  reaction  chamber. 
Distilled  water  is  then  allowed  to  enter  the  reaction 
chamber.  Acetylene  is  rapidly  produced  and  captured  as 
a  white  solid  by  trap  C. 

The  Eigh-Frequency  Cracking  of  Acetylene 

The  production  of  carbon  targets  by  a  high-frequency 
(h,  f . )  discharge  is  described  in  the  references  (53), 
The  apparatus  used  in  the  present  experiment  was  designed 
so   that  it  would  rest  on  the  base  of  a  conventional 
evaporation  system  (thereby  using  the  diffusion  pump,  cold 
trap  and  associated  electronics  of  the  evaporator)  (figure  13), 
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Figure  13.   High-frequency  discharge  apparatus. 
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The  target  backings  (10  mil  Ta)  were  placed  in  the  3^.5  cm 
long  by  3.8  cm  I.  D.  "Pyrex"  tube.   The  tube  could  accommodate 
four  backings  (three  resting  in  the  horizontal  tube  and  one 
attached  to  the  endplate)  without  affecting  the  cracking 
process.   A  Tesla  coll  furnished  the  h.  f.  discharge,   The 
system  is  evacuated  to  a  pressure  of  about  3  x  10""''  torr, 
then  isolated  from  the  diffusion  pump.   The  Tesla  coil  is 
activated,  and  acetylene  is  admitted  into  the  cracking 
chamber  through  the  "side-arm"  from  trap  C  (which  was  removed 
from  the  apparatus  used  to  produce  the  acetylene).   The  C2H2 
was  allowed  to  sublime  by  removing  the  trap  from  the  liquid 
nitrogen  and  placing  it  in  an  alcohol-dry  ice  slush,  thereby 
allowing  acetylene  to  enter  the  cracking  chamber.   To  make 
uniform  targets  required  careful  control  of  the  C2H2  vapor 
pressure.   This  was  done  by  using  a  stopcock  to  limit  the 
flow  of  C2H2.  Experience  showed  that  the  most  efficient 
pressure  could  be  determined  by  the  color  of  the  glow  dis- 
charge.  This  optimum  pressure  is  indicated  by  a  pale  green 
glow  with  no  visible  striations.  After  the  acetylene  has 
completely  cracked,  the  color  of  the  glow  discharge 
changes  to  blue,  characteristic  of  nitrogen.  At  this  point 
the  process  could  be  halted  (producing  very  thin  targets) 
or  repeated  several  times  (producing  thicker  targets).   It 
was  observed  that  the  target  attached  to  the  endplate  was 
thicker  than  the  ones  resting  in  the  glass  tube.   It  was 
found  that  by   starting  with  20  mg  of  enriched  12c  it  was 
possible  to  produce  eight  targets.   Four  were  thin,  three 
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were  of  intermediate  thickness  (15-30  keV  for  3.6  MeV  alpha 
particles),  and  one  was  relatively  thick  (about  55   keV  for 
3.6  KeV  alpha  particles). 

Targets  prepared  by  this  procedure  were  able  to  with- 
stand beam  currents  of  about  3  /xA;  furthermore,  the  spectra 
Indicated  a  significant  enrichment  in  12C  (see  next  section). 
The  advantage  of  this  procedure  was  that  after  the  method 
was  perfected  a  sizable  number  of  enriched  12C  targets 
was  available  at  minimal  cost.   Of  course,  the  above  procedure 
could  also  by  employed  to  produce  targets  enriched  in  isotopes 
of  carbon  other  than   C. 


CHAPTER  III 


EXPERIMENTAL  RESULTS 


The  yield  curves  which  follow  were  obtained  with 
the  12.7  era  x  12.7  cm  Nal(Tl)  crystal  fixed  at  a  S5°   angle 
with  respect  to  the  alpha  particle  beam.   This  detector  was 
usually  0.5  cm  from  the  target;  however,  some  long  runs  were 
obtained  with  the  detector  at  2.5  cm  from  the  target.   Since 
the  detector  was  quite  close  to  the  target,  summing  of 
cascade  gamma  rays  may  occur.  A  20.1  cm3  Ge(Li)  detector 
was  also  employed  in  examining  the  7  MeV  cascade  radiation 
from  the  9.85  MeV  level  in  l60. 

The  target  thicknesses  were  determined  by  employing 
the  narrow  level  (Tcm  =  0.75  keV)  in  l60  at  an  excitation 
of  9.85  MeV  (26,  27).   This  part  of  the  excitation  function 
for  the  target  used  for  the  9.85  MeV  measurements  is  shown 
in  figure  14.   No  corrections  have  "been   made  to  this  curve. 
By  determining  the  full  width  at  half -maximum,  the  target 
thickness  was  found  to  be  25  keV  for  3.6  HeV  alpha  particles 
incident  at  55°   to  the  surface.   Using  standard  methods  the 
thicknesses  were  extrapolated  to  the  lower  alpha  particle 
energies  as  required.  A  thicker  target  {$$   keV  for 
3.6  MeV  alpha  particles  incident  at  $5°)   was  used  for  most 
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Figure  14.  Enriched  12C  excitation  function  in  the 
region  of  the  Ea  =  3.58  MeV  resonance. 
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off -resonance  cross  section  measurements.   Some  off -resonance 
cross  sections  were  also  determined  with  a  30  keV  (for 
3.6  KeV  alpha  particles)  target. 


Exci  tat  ion  Fung  t  i  ons 


The  excitation  functions  in  the  region  E  =  2.7-3.7  KeV 
for  targets  of  enriched  12C ,  natural  carbon  and  enriched  13c 
are  shown  in  figures  15,  16  and  17,  respectively.   The 
ordinate  in  figures  15  and  16  represents  the  number  of  counts 
per  300  n  Coulombs  of  singly  charged  helium  ions  within  the 
region  Ey  =  8.7  MeV  to  10.8  KeV.   In  the  excitation  function 
for  the  enriched  13c  target  the  ordinate  represents  the 
number  of  counts  (in  the  same  sum  region  as  for  the  natural 
carbon  and  enriched  12C  targets)  per  150 n  Coulombs.   Beam 
currents  were  usually  between  1.5-3  nA.   The  detector  was  at 
a  distance  0.5  cm  from  the  target  and  at  an  angle  of  55° 
with  the  beam  axis.   The  enriched  12C  target  was  15  pg/cm2 
thick.,  the  natural  carbon  target  was  22  jzg/cm2  thick  and 
the  enriched  (probably  l5-25  per  cent  enrichment)  13C  target 
was  Ik  Mg/cm2  thick.   The  enrichment  of  *2C  is  indicated  by 
comparing  the  resonance  at  Eflt  =  3.58  KeV  in  the  three  figures. 
The  broad  peak  in  the  enriched  12C  spectrum  at  Ea  =  3.24  KeV 
is  in  the  region  of  the  9.59  KeV  excitation  energy  in  160| 
however,  most  of  this  yield  is  probably  due  to  the  ^cCa.n^O 
reaction.   It  is  difficult  to  determine  the  actual  enrichment 
of  12C,  since  some  of  the  background  is  due  to  neutrons 
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originating  at  collimators,  etc.   Certain  runs  on  tantalum 
backings  (at  Ea  =  3,  24  MeV)  indicate  that  this  source  of 
background  may  be  as  much  as  30  per  cent  of  the  total  back- 
ground. 

An  indication  of  a  resonance  due  to  the  12C(a,Y)l60 
reaction  at  Ea  =  3. 48  MeV  was  observed  in  the  excitation 
functions  for  the  enriched  12C  and  natural  carbon  targets. 
This  structure  was  not  observed  in  the  excitation  function 
for  the  enriched  Y\  target.  No  attempt  was  made  to  determine 
conclusively  whether  this  structure  could  be  attributed  to 
the  12C(a,Y)l60  reaction. 

The  9.86  MeV  Lege]  In  l6" 

The  level  at  Ea  =  3.58  MeV  corresponds  to  a  level  in 
1  0  at  9.85  MeV  (2+).   This  level  has  been  observed  through 
the  12C(a,Y)l60  reaction  by  Meads  and  Mclldowie,  and  Larson  and 
Spear  (10,  12).   A  spectrum  taken  at  Ea  =  3.58  MeV  is  shown 
in  figure  18.   No  corrections  have  been  made.   Larson  had 
observed  7   MeV  cascade  radiation  in  the  spectrum  from  his 
thick  target  (11).   This  7  MeV  radiation  is  also  visible  in 
figure  18.  Because  of  the  resolution  of  the  Nal(Tl)  detector 
it  is  difficult  to  determine  the  level  in  l60  through  which 
this  cascade  proceeds.   For  this  reason  the  level  was 
examined  using  a  20.1  cm3  Ge(Li)  detector.   A  spectrum  from 
the  1024  channels  of  the  TMC  analyzer  is  shown  in  figure  19. 
This  is  a  point-by-point  sum  of  individual  runs  of 
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Figure  18.   Gamma  ray  spectrum  from  Nal(Tl)  detector  at 

Ea  =  3.58  MeV.   The  three  arrows  correspond  to  the  positions 

of  the  photopeak  and  the  two  escape  peaks. 
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0,09,  0.05  and  0.12  Coulomb.   The  (9.85  MeV-2e)  peak  is 

clearly  visible.   The  smaller  peak  is  from  the  (6.92  MeV-2e) 

radiation.   The  peak  corresponding  to  the  2+ 2+  transition 

((9.35  -  6.92)  MeV-2e)  is  also  evident.   The  other  peaks 
apparently  correspond  to  various  background  gamma  rays.   From 
this  it  appears  reasonable  that  the  cascade  radiation  proceeds 
through  the  2+  level  at  6.92  MeV  in  l60.   The  results  for 
the  radiative  widths  (determined  from  Nal(Tl)  data)  for  the 
9.85  MeV  ground  state  and  cascade  transitions  are  given  in 

table  3.   The  value  T     =  0.0050  eV  corresponds  to 

2  9  o 

IMI   =  2.76  x  10"*,  where  |Hl  z  is  the  ratio  of  V     to  the 

radiative  width  calculated  from  an  extreme  single  particle 
model  (Tw).  Tw  is  the  so-called  Welsskopf  limit.   In  deter- 
mining the  radiative  x*idth  for  the  cascade  radiation,  Larson 
assumed  an  isotropic  distribution  since  he  was  not  able  to 
determine  the  identity  of  the  intermediate  level  (11). 

The  sources  of  error  are  given  in  table  k.      Errors 
due  to  current  integration,  extrapolation  of  stopping  cross 
sections  to  different  energies,  and  geometry  have  been 
combined  and  listed  in  the  column  labeled  "other".   The  total 
error  is  the  square  root  of  the  sum  of  the  squares  of  the 
individual  errors.   From  table  k   it  can  be  seen  that  no  error 
is  listed  under  "target  thickness"  for  the  9.85  MeV  radiations. 
Due  to  the  extremely  narrow  total  width  of  this  state 
(rcm  =  °«75  keV)  the  factor  arctan(6/T)  in  the  expression 
for  the  radiative  width  reduces  to  tt/2  and  is  not  a  source 
of  error.   The  difficulty  in  extracting  the  9.85  MeV  cascade 
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radiation  In  the  presence  of  the  strong  ground  state 

radiation  caused  a  10  per  cent  error  for  "background  sub- 
traction" for  the  cascade  radiation.  The  brancing  ratios 
for  the  9.85  MeV  level  are  given  in  figure  20. 

Cross  Sections  and  the  9.59  MeV  Level 

The  1"  state  at  9.59  MeV  excitation  in   0  was 
exaimined  with  the  29  /ig/cm2  thick  target.   This  level  la 
very  broad,  having  a  center-of-mass  width  of  6^5  keV  (26). 

Figure  21  shows  a  spectrum  from  the  Kal(Tl)  crystal  at 
E  =3.2^  MeV.   There  were  no  corrections  made  to  this 
spectrum.   The  detector  was  0.5  cm  from  the  target  and 
at  an  angle  of  55°.      The  total  amount  of  charge  collected 
was  5.5 x  10  n Coulombs  of  singly  Ionized  helium.   The  arrows 
indicate  the  expected  positions  of  the  photopeak  and  escape 
peaks.   For  comparison  purposes  a  spectrum  at  Ea  =  2.95  MeV 
obtained  from  the  enriched  -*C  target  (Ik   fxg/cm2)  is  shown 
in  figure  22.   Spectra  obtained  at  E  =  3.2^  MeV  with  this 
target  have  the  same  general  appearance.   The  spectrum 
from  the  ^C   is  exponentially  decreasing  under  the  radiation 
from  the  9.59  MeV  level.   This  same  feature  is  also  observed 
In  the  spectra  from  a  natural  carbon  target  (figure  23). 

The  9.59  MeV  yield  was  obtained  by  subtracting 
fractional  parts  of  the  9.85  MeV  spectrum  from  the  9.59  MeV 

spectrum.   The  resulting  spectrum  was  then  compared  with  a 

13      1 6> 

C(a»n)   0  spectrum  to  determine  its  fit  with  the  background. 
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Figure  20.   Branching  ratios  of  9.85  MeV  level  in  l60. 
The  energy  levels  are  not  drawn  to  scale. 
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Figure  21.  Gamma  ray  spectrum  from  Nal(Tl)  detector  at 
Ea  =  3.24  MeV  for  enriched  12C  target. 
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Figure  22.  Gamma  ray  spectrum  at  Ea  =  2.95  MeV  for  enriched 
1JC   target. 
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Figure  23.   Gamma  ray  spectrum  at  E  =  2.95  MeV  for  natural 
carbon  target. 
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This  same  method  was  also  used  in  determining  the  off- 

resonance  yields.   This  procedure  gave  a  value  of 

0,021  +  0.004  eV  for  the  radiative  width.  A  comparison  with 

previous  results  is  given  in  table  5.   The  present  work 

agrees  well  with  the  value  given  by  Larson  (11).   The  value 

TY  =  0.021  eV  gives  a  value  of  5.5   x  10"-5  for  |H|2. 

Although  there  is  some  indication  of  7  MeV  cascade 

radiation  in  the   C  spectrum  near  channel  280,  this  feature 

is  also  seen  in  the  ^C  spectrum.   It  is  felt  that  the 

spectrum  peak  in  channel  280  is  the  result  of  the  capture 

of  thermal  neutrons  by  the  *%a  in  the  Nal(Tl)  crystal. 

When  a  paraffin  shield  was  Installed  between  the  detector 

and  the  collimators*  the  spectrum  peak  became  more  pronounced. 

The  curves  shown  in  figures  21,  22  and  25  were  obtained 

with  this  shield,  while  the  data  shown  in  figures  23  and  2k 

were  obtained  without  the  paraffin  shield.   There  was  some 

indication  that  although  this  shield  increased  the  7  MeV 

background,  it  decreased  the  background  in  the  9  MeV  region. 

It  has  been  observed  that  9  MeV  radiation  (resulting 

12      1^ 
from  the   C(d,y)  N  reaction)  may  also  present  a  background 

problem  in  a   C  +  a  experiment  (11).   The  reason  presented 

is  that  if  deuterons  exist  in  the  ion  beam,  the  beam  energy 

analyzing  magnet  will  permit  the  target  to  be  bombarded  with 

deuterons  having  one-half  the  energy  of  the  alpha  particles, 

since  molecular  D£  ions  will  have  nearly  the  same  path  as 

h  12 

singly  ionized  He.   The  bombardment  of   C  by  1.65  MeV 

deuterons  produces  9.^  MeV  cascade  radiation  to  the  2.31  MeV 
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81 
first  excited  state  In   N.   The  r.  f.  source  bottle  used 

for  this  experiment  has  never  been  used  to  Ionize  any 

hydrogen  isotope?  hence,  it  was  assumed  that  the  contribution 

from  deuteron-induced  reactions  was  negligible  in  the  present 

experiment. 

The  off -resonance  cross  sections  were  determined  from 

spectra  such  as  the  one  shown  in  figure  2k .   The  52  MS/cm6 

thick  target  was  used  for  this  data.   The  ground  state 

radiation  is  clearly  visible  in  channels  355-390.  A 

comparison  with  figure  22  shows  that  the  number  of  counts  in 

13 

the  corresponding  channels  for  the  ^C  spectrum  is  decreasing 

almost  exponentially.  For  the  spectrum  at  Ea  =  2.95  MeV 
(as  for  the  other  cross  section  measurements)  the  12,7  cm  by 
12.7  cm  Nal(Tl)  detector  was  at  a  55°   angle  with  respect  to 
the  beam  axis,  and  i*as  located  0.5  cm  from  the  target.   The 
ground  state  radiation  is  more  visible  at  E  =  2.95  MeV 
(figure  2k)    than  at  the  peak  Ea  =  3.2^  HeV  (figure  21).   Two 
factors  contribute  to  this  appearance.  First,  the  average 
number  of  counts  per  channel  in  figure  21  Is  about  10  times 
that  In  figure  23;  hence,  on  the  logarithmic  scale  the  peak 
is  not  as  evident.  Also,  the  ^C(atn)   q  cross  section  is 
about  5  times  as  large  at  Ea  =  3.2^  HeV  as  at  E  =2.95  MeV 
(kj) ,      Thus,  the  ground  state  radiation  is  more  difficult 
to  observe  at  the  resonant  energy  because  of  the  increased 
neutron  background. 

The  yield  spectra  for  the  other  data  points  are 
similar  to  the  spectrum  shown  in  figure  2k.     As  the  alpha 
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particle  bombarding  energy  decreases,  the  spectra  resemble 
the  -^C  background  more  closely.   For  example,  at  Ea  =  2.1^  MeV 
there  Is  only  a  slight  indlcabion  of  the  ground  state  radiation. 
This  is  shown  in  figure  25 »  where  the  data  were  obtained  with 
the  28.8  pig/ cm2  thick  target.   Because  of  this  the  error 
caused  by  extracting  the  yields  from  the  backgrounds  are 
relatively  large  for  the  cross  section  meaurements  made  at 
low  bombarding  energies  (table  *0.   The  amount  of  singly 
ionized  helium  collected  was  typically  between  0,0^5  and 
0.16  Coulomb. 

The  yield  at  each  of  the  experimental  points  was 
normalized  to  that  at  3.58  MeV  by  multiplying  the  total  number 
of  counts  in  the  capture  spectrum  peak  by  the  ratio  of  the 
number  of  channels  enclosing  the  peak  to  the  number  of  channels 
in  the  standard  spectrum  peak.   The  values  of  the  cross 
section  are  listed  in  table  6.   The  errors  shown  were  computed 
from  the  values  given  in  table  k. 
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Table  6 
The  C(a,y)     0  reaction  cross  section. 


He  energy  Cross  section 

(MeV)  (10-9  barns) 


3.37  32.0  +  6.1 

3.2^  35.0  +  6.1 

3.22  35.0  +  6.1 

3.10  29.0  +  5.8 

3.07  19.8  +  3.8 

2.95  15.2  +  3.0 

2.92  11.9  ±  2.4 

2.76  8.0  +  1.8 

2.5^  4.3  ±   1.3 

2.35  3.5  +  1.2 

2.14  1.5  +  0.5 

1.90  0.3  +  0.1 


CHAPTER  IV 


DISCUSSION 


'Resonances 


In  this  experiment  the  radiative  capture  of  alpha 
particles  by  12C  below  Ea  =  3.6  MeV  has  been  examined.   The 
7  MeV  cascade  radiation  from  the  9.85  MeV  (2+)  level  in   0 
was  determined  to  proceed  through  the  6.92  MeV  (2+)  level. 

It  would  be  of  interest  to  examine  this  (0+ ^  2*  M1E2  >  2+) 

transition  using  angular  distribution  methods  to  determine 
the  mixing  parameter  X,  where  X  a  S(M1)/S(E2).   S(M1)  and 
S(E2)  are  amplitudes  for  dlpole  and  quadrupole  emission, 
respectively  (29).   Further  study  of  this  state  would 
probably  best  be  accomplished  through  angular  correlation 
methods.  At  the  time  this  experiment  was  performed,  a  new 
correlation  table  was  being  designed  and  constructed  by 
members  of  the  laboratory.   This  table  will  be  able  to  make 
the  above  measurements  accurately.   For  this  reason  no 
determination  of  the  mixing  parameter  X  was  attempted  using 
the  present  apparatus. 

The  radiative  width  measurements  agree  in  general 
with  the  results  of  Larson  (the  errors  overlap)}  however, 
the  values  tended  to  be  slightly  smaller  than  the  values 

8^ 
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obtained  by  Larson.   This  may  be  due  in  part  to  the  method 

used  In  determining  the  peak-to-total  ratio  for  the  present 

experiment. 

Target  Deterioration 

Target  deterioration  was  present  only  at  higher  beam 
currents  (about  k~7  pA-)>'    At  these  beam  currents  the  carbon 
was  dissipated,  and  the  tantalum  backing  became  quite  pitted. 
Carbon  deposition  was  not  observed  on  the  target,  even  after 
the  relatively  long  bombarding  times.   In  order  to  observe 
this  more  carefully,  several  relatively  long  runs  (greater 
than  0.04  Coulomb)  were  made  on  clean  Ta  backings  using 
both  low  (1-3  MA)  and  high  (4—8  jzA)  beam  currents.   No 
deposition  of  organic  material  was  observed  for  either  of 
these  currents.  However,  the  higher  beam  currents  did  cause 
the  tantalum  to  change  color  at  the  spot  which  had  been 
bombarded  by  the  beam.   This  spot  was  gold  in  color  on  both 
sides  of  the  backing.  This  is  believed  to  be  caused  by  the 
tempering  of  the  tantalum  by  the  concentration  of  heat  on  the 
beam  spot.  This  same  effect  was  observed  during  preliminary 
experiments  in  target  preparation  by  thermal  cracking  methods. 
When  the  Ta  backing  was  outgassed  by  heating  the  backing  to 
a  bright  yellow  glow  in  a  vacuum,  the  Ta  backing  would  have 
the  characteristic  gold  color  and  be  very  rigid  when  removed} 
hence,  it  seemed  that  the  backing  had  been  tempered  by  the 
high  temperature  used  in  the  outgassing  process.   It  might 
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be  noted  here  that  there  was  some  evidence  (using  the  natural 
carbon  target  made  by  h.  f .  cracking)  to  Indicate  that 
backings  that  had  been  outgassed  in  a  vacuum  produced  targets 
which  exhibited  greater  resistance  to  deterioration  (by 
high  beam  currents)  than  targets  produced  from  backings  that 
had  not  been  outgassed.  Enriched   C  targets  were  made  using 
Cu  backings,  but  these  were  not  used  in  this  experiment. 
These  may  be  able  to  withstand  higher  beam  currents  than 
those  made  with  tantalum  backings. 

Cross  Sections 

The  cross  section  for  the  12C(a,Y)   0  reaction  at  low 
alpha  particle  bombarding  energies  is  necessary  for  an  accurate 
extrapolation  to  effective  stellar  energies.   The  Coulomb 
barrier  B  for  the  12C  +  He  can  be  approximated  by 


Bc  =  ZiV2/5' 


—     —         to 

where  R  =  1.4  x  (Ax3  +  A03  )  x  10"13  cm. 

The  prescription  gives  the  value  3.18  MeV  for  B  ;  hence, 

below  3  MeV  barrier  penetration  becomes  increasingly  Important. 

Tombrello  has  computed  the  effect  of  Interference 
between  the  1"  states  at  7.12  MeV  and  9.59  MeV  in  1^0. 
Using  a  value  ea  =0.10  for  the  reduced  width  of  the  7.12  MeV 
state,  he  has  calculated  the  values  of  the  capture  cross 


8? 
section  for  both  constructive  and  destructive  interference 
(22).  Since  the  signs  of  the  amplitudes  were  unknown,  it 
was  Impossible  to  determine  whether  the  interference  was 
constructive  or  destructive. 

The  values  for  the  12C(a,v)   0  reaction  cross  section 
obtained  from  the  present  experiment,  and  those  computed  by 
Tombrello  (assuming  ea2  =0.1)  are  plotted  in  figure  26. 
The  value  obtained  for  the  capture  cross  section  at 
Ea  s  3  KeV  is  slightly  smaller  than  that  calculated  by 
Tombrello;  while  at  Ea  =  1.9  MeV  it  Is  about  1.5  times  as 
large  as  the  calculated  value.   The  results  of  this  experi- 
ment indicate  that  the  interference  is  constructive  and  that 

p 

0a  is  between  0.1  and  O.J, 


The  12c(atV)l60  Reaction  Rate 
at  Effective  Stellar  Energies 


The  reduced  alpha  particle  width  ea2  of  the  7.12  MeV 
level  in  160  is  essential  in  the  determination  of  the 
12C(a,Y)   0  reaction  rate  at  stellar  temperatures  (see  the 
section  on  theory).   This  reaction  rate  is  important  in 
stellar  model  calculations  since  it  determines  the  relative 
abundances  of  X*C  and   0  and  hence  the  futher  evolution  of 
the  star. 

This  experiment  indicates  the  cross  section  factor  S0 
at  a  typical  stellar  temperature  of  108  °K  has  the  value 

So(12C(a,y)l60)  =  (0.65  ±  0,32;  x  103  keV  barns, 
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Figure  26.   The   c(a,Y)   0  reaction  cross  section  values. 
Ihe  solid  line  has  been  drawn  through  the  calculated  points 
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2 
using  9   =0.2+0.1.   Slmllary  (using  the  equations  derived 
a       ~ 

in  the  section  on  theory),  the  present  experiment  yields 
the  following  value  for  the  12C(a,y)   °  reaction  rate  at  a 

a 

stellar  temperature  of  10  °K: 


p  (12C)  =  (0.84-  +  0.^2)  x  10~20X  Pf._   s"1, 

where  P  is  the  density  of  the  core  in  g  cm"*3  (~10^) ,  X  is 

.     "  a 

the  abundance  (~1)  by  weight  of  He  and  f19  is  the  screening 

factor  (37).  Although  the  above  errors  on  S  and  p  may 

appear  large,  it  should  be  noted  that  prior  to  this  experiment, 

2 

9   could  only  be  estimated  to  within  a  factor  of  10  at  best. 

Also,  it  is  quite  probable  that  the  errors  will  be  further 
reduced  as  revised  fits  are  obtained  to  the  measured  values 
of  1;he  capture  cross  section. 

A  determination  of  12C  +  He  capture  cross  sections 
at  lower  bombarding  energies  than  those  employed  in  this 
experiment  would  also  produce  a  better  value  for  the  reaction 
rate.  Experimental  procedures  that  could  reduce  the 

-^C(a,n)   0  background  would  be  necessary.   Using- a  more 
highly  enriched  12C  target,  and  time-of -flight  methods  to 
eliminate  neutron  background  originating  from  the  collimators 
may  be  the  most  useful  procedure  since  it  appeared  that  some 
of  the  neutron  background  resulted  from  the  beam  striking 
collimators  and  slits  on  which  carbon  had  deposited. 


APPENDIX 


GAMKA  HAY  DETECTION  EFFICIENCIES  OF  Nal(Tl)  CRYSTALS 

Total  Incident  Intrinsic  Sfflcleny 

The  "total  incident"  intrinsic  efficiency,  l^(Ey)t    of 
a  Nal(Tl)  crystal  in  the  form  of  a  right  cylinder  for  gamma 
rays  of  energy  E  is  defined  as  the  probability  for  the 
detection  of  a  gamma  ray  that  is  incident  on  the  face  of  the 
crystal  (56>   57).      This  is  equal  to  the  ratio  of  the  total 
number  of  gamma  rays  detected  per  second  to  the  total  number 
of  gamma  rays  incident  on  the  face  of  the  crystal.   In  the 
form  of  an  equation,  this  Is 

NT 
*T(E)  = 


N0f 

where  N<p  is  the  total  number  of  gamma  rays  of  energy  E 
detected  per  second,  N0  is  the  total  number  of  gamma  rays 
of  energy  E  emitted  per  second  by  an  Isotropic  source,  and 
f  is  the  fraction  of  kv   that  the  solid  angle  of  the  crystal 
face  subtends.   Stated  differently,  f  is  defined  by 


r.-i 


where  ft  is  the  solid  angle  subtended  by  the  crystal. 
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It  Is  evident  that  %/EL  Is  simply  the  fraction  of  the 
gamma  rays  detected,  and  f  is  simply  the  probability  that  a 
gamma  ray  from  an  isotropic  source  will  strike  somewhere  on 
the  face  of  the  crystal.  The  geometry  for  the  case  of  an 
isotropic  point  source  of  gamma  rays  of  energy  E  located 
on  the  axis  of  a  right  cylinder  of  radius  R  and  thickness  t 
is  shown  in  figure  27. 

The  parameter  x2(e')  is  the  path  length  in  the  Nal(Tl) 
crystal.   The  parameter  x-^e*)  is  the  path  length  in  any 
absorbing  material  preceding  the  crystal.   The  quantities 
x?  and  Xn  have  the  following  explicit  dependency  on  e'i 


x^( 6 ' )  =  d  sece ' ,  for  0  <  e '  <  6  =  arctan 


d 

R 


x«(e')   =  t  sece',   for  0<  e'  <  e,    =  arctan 


t  +  d 


xo(e')   =  R   csca'    -  d  sece',    for   ei <  e'<  ert  =  arctan  . 

1      °         d 


The  intensity  of  the  source  will  be  attenuated  by  a 
factor  exp(-  M^(E)x-.  (e  ' ) )  in  traversing  the  absorbing  material 
preceding  the  crystal.   (i,  is  the  total  linear  absorbing 
coefficient  for  a  gamma  ray  of  energy  E  in  the  absorbing 
material  (58).   In  traversing  the  Nal(Tl)  crystal  the  gamma 
ray  intensity  will  be  attenuated  by  a  factor  exp(-  M 2(E)x2( 9  * )) » 
where  ^2(E)  *s  the  total  linear  absorption  coefficient  for 
a  gamma  ray  of  energy  E  in  Kal(Tl). 

N  /4tt  is  the  intensity  per  unit  solid  angle  (the  number 
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a 


o 


■p 
as 


o 

H 


■P 


9^ 

of  disintegrations  per  second  per  unit  solid  angle)  of  the 
isotropic  source.   The  number  of  gamma  rays  emitted  per 
second  into  the  small  solid  angle  AQ'  about  £2'  is  (N0/4rr) A  Q' . 
In  the  small  solid  angle  Afi' the  number  of  gamma  rays  traversing 
the  crystal  without  interacting  is  equal  to  the  quantity 


ANn  =   f_  exp(-M1(E)x1(e')  -  M2(E)x2(9'))ASr. 


^rr 


Hence,  the  total  number  N  of  gamma  rays,  emitted  by 
an  isotropic  source  of  intensity  N  not  interacting  with 
the  crystal  is  given  by 

^-tt  ^crystal 

The  number  Ninc  of  gamma  rays  incident  on  the  face  of  the 
crystal  is  given  by 

^tt    crystal 

Therefore,  NT  (the  number  of  gamma  rays  detected  per  second, 
assuming  that  each  interaction  produces  an  observable  pulse) 
is  obtained  from  the  following: 

NT  =  Ninc  -  Nn»  or 

T  ~  -^  Je  (1  -  e    ^    2     )  d0.t 
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Using  the  relation  derived  for  NT,  the  total  incident  intrinsic 

efficiency  T r  is  given  by 

„  (EJ  «JL  fNl  -  e-^^je-^^l  sine'de', 
1  Y    2f  JQ 

where  f  =  fi/4n  =  (IAtt)  f  d.Q'   =  |(l-d(R2  +  d2)"2), 

•'crystal 

It  is  often  convenient  to  define  a  "source"  intrinsic 
efficiency,  T?fc(E),  as  the  product  of  f  and  Tr 
(>7t  =  fi*T  =  NT/N0).   In  this  case  the  geometry  factor  f  is 
accounted  for  in  the  computation  of  the  source  intrinsic 
efficiency.  Hence,  the  absolute  intensity  NQ  of  a  source 
is  obtained  by  dividing  the  observed  counting  rate  N.p  by 
i?t.   If  we  introduce  Legendre  polynomials  Pn(cose')»  we  will 
be*  able  to  relate  the  source  intrinsic  efficiency  17^  to  the 
integrals 

Jn  =  __!_  f   Pn(cose')e"  "1*1(1  -  e"  M2X2)  dfi'. 
2rr  "'crystal 

These  are  the  integrals  defined  by  Rose  (30).   Therefore, 

NT     JQ 
i|t(E)  = - 

No     2 

is  the  source  intrinsic  efficiency. 

In  general  M2  is  much  larger  than  f*^,  and  it  is  possible 
to  calculate  values  of  17  ^  using  ft^  -   0«   Values  of  1^ 
(computed  at  Oak  Ridge  with  ^  *  0)  and  values  of  *Jt 
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(computed  by  Wolickl  with  /x,  =  0)  are  compiled  in  Kott  and 
Sutton  for  several  different  crystal  sizes  and  source- 
detector  geometries  (56,  59).  For  other  geometries  and 
crystal  sizes  one  may  interpolate  between  given  values.   More 
curves  for  ^  may  be  found  in  the  article  by  Rivers  (60). 
These  calculations  have  been  further  extended  by  Heath  to 
include  disc  and  line  sources  (6lf  62).   Values  of  absorption 
coefficients  for  photons  have  been  compiled  by  Grodstein  (63). 

The  pulse  spectrum  for  monoenergetlc  gamma  rays  is 
complicated  by  pulses  from  externally  scattered  photons  in 
addition  to  the  spectrum  resulting  from  the  primary  photons 
of  the  source.   The  scattered  radiation  results  from  the 
Compton  gamma  rays  scattered  into  the  crystal  by  shielding, 
wall,  target  or  source  holders,  etc.   This  radiation  has 
a  peak  which  occurs  at  an  energy  corresponding  to  180  Compton 
scattering,  and  hence  is  due  mainly  to  backscattered  photons 
emerging  from  the  photomultiplier  window  (5&t   58). 

Thus,  in  order  to  measure  absolute  yields  the  spectrum 
of  pulses  must  be  corrected  accordingly.   In  general  it  is 
relatively  easy  to  obtain  the  yield  for  monoenergetlc  gamma 
rays.   However,  it  is  much  more  difficult  to  accurately 
determine  the  absolute  yields  if  two  or  more  gamma  rays 
of  different  energies  are  present.   The  analysis  of  a  spectrum 
containing  gamma  rays  of  several  energies  is  accomplished 
by  "spectral  decomposition"  (6k).      That  is,  the  spectrum  is 
broken  into  its  individual  components.   After  corrections 
have  been  made  for  background  and  backscattering,  the 
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"response  function"  (the  detailed  shape  of  the  spectrum  of 

a  single  gamma  ray  energy)  of  the  gamma  ray  having  the 
largest  energy  is  obtained  and  subtracted  from  the  spectrum. 
The  process  is  repeated  for  the  residual  spectrum.   This  is 
repeated  until  the  response  functions  for  all  of  the  gamma 
rays  have  been  obtained.   These  individual  "component  spectra" 
may  then  be  analyzed  separately. 

Source  Full  Energy  Peak  Efficiency 

Since  a  pulse  height  spectrum  depends  so  strongly  on 
the  particular  experimental  arrangement,  it  is  desirable  to 
choose  a  portion  of  the  spectrum  that  is  somewhat  less 
sensitive  to  scattered  radiation  and  background  (56,   59,  6k) . 
Now,  the  shape  of  the  photopeak  is  relatively  insensitive 
to  radiation  scattered  into  the  Nal(Tl)  crystal,  and  back- 
ground corrections  are  usually  small.   Furthermore,  the 
shape  of  the  photopeak  is  very  similar  to  a  normal  curve 
(Gaussian) ,  and  thus  the  total  number  of  counts  in  the 
peak  (which  is  equal  to  the  area  under  the  peak)  may  be 
related  to  certain  parameters  of  the  Gaussian  (peak  counting 
rate,  full  width  at  half- maximum,  energy  Interval  per 
channel)  (65,    66).      Hence,  the  problems  encountered  in 
relating  an  observed  counting  rate  to  the  absolute  intensity 
N0  of  a  given  source  are  greatly  simplified  if  one  employs 
the  total  number  Np  of  counts  in  the  full  energy  peak  instead 
of  employing  Nm, 
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The  "incident"  intrinsic  peak  efficiency  (not  the 
source  full  energy  peak  efficiency  17(E)),  may  now  be  defined 
as  the  probability  that  a  gamma  ray  of  energy  E  incident  on 
the  Nal(Tl)  crystal  will  produce  a  pulse  lying  in  the  full 
energy  peak.   Stated  differently,  the  "incident"  intrinsic 
peak  efficiency  is  the  ratio  of  the  number  of  counts  per 
second  lying  in  the  photopeak  to  the  number  of  gamma  rays 
of  energy  E  incident  per  second  on  the  face  of  the  crystal 
(assuming  i^1   =  0).   Thus,  the  "incident"  intrinsic  peak 
efficiency  is  given  by 

"peak  -  VNof ' 

It  is  now  possible  to  define  the  source  full  energy 
peak  efficiency  77p  as  a  product  of  f  and  V-peak* 

9p(E)  =  f^^   =  Np/N0. 

It  is  quite  difficult  to  calculate  rj     directly  from 
the  theory  because  a  pulse  lying  in  the  photopeak  may  have 
been  caused  by  the  summing  of  several  successive  processes 
in  the  stopping  of  a  single  gamma  ray.   For  example,  a 
photon  may  be  absorbed  by  a  Compton  scattering  followed  by 
pair-production  followed  by  annihilation  of  the  positron. 
If  all  quanta  and  electrons  are  absorbed,  the  resulting 
light  from  the  scintillator  will  produce  a  pulse  lying  in 
the  photopeak.   This  is  possible  since  the  time  required  for 
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collecting  the  light  is  long  compared  to  the  time  required 

by  the  scintillator  to  absorb  the  gamma  ray.   This  pulse 
would  then  be  indistinguishable  from  one  due  to  a  single 
photoelectric  process.   Hence,  it  is  extremely  difficult 
to  determine  the  actual  absorption  coefficient  /x  that  would 
be  required  in  the  calculation  of  tj  .   However,  by  defining 
a  suitable  ratio  R  it  is  possible  to  relate  r;  to  the 
calculated  source  intrinsic  efficiency  17^.   The  ratio  R  is 
determined  experimentally  as  a  function  of  the  gamma  ray 
energy  E  .  R(E)  is  defined  as  the  ratio  of  the  number  Np 
of  counts  under  the  photopeak  (area  of  photopeak)  to  the 
total  number  N™  under  the  entire  pulse  spectrum  (corrected 
for  background,  backseat tering,  and  spectrally  decomposed). 
Hence,  the  relation  between  the  source  full  energy  peak 
efficiency  T?p  and  *?£  is  given  by 


np  =  R  V 


The  "peak-to-total"  ratio  R  may  be  determined  by  analyz- 
ing spectra  from  standard  sources  (if  the  energy  of  interest 
is  low).   If  the  energy  range  of  the  gamma  rays  is  higher, 
then  one  must  resort  to  certain  nuclear  reactions  where 
there  is  a  predominant  gamma  ray  present  in  this  range.  In 
this  experiment  the  9.17  MeV  gamma  ray  resulting  from  the 
13c(p,y)   N  reaction  was  employed  (figure  28).   The  9.85  MeV 
gamma  ray  resulting  from  the  12C(cx,y)   0  reaction  was  also 
used  in  determining  R, 
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Figure  28.   Gamma  ray  spectrum  from  the  13c(p,Y)  ^  reaction 
(Ey  =  9.17  MeV).   The  dashed  line  indicates  the  horizontal 
extrapolation  used  in  the  present  experiment  to  determine 
the  peak-to-total  ratio  R.   This  gives  the  value  E  =  0  50 
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Although  the  spectrum  exhibits  a  strong  gamma  ray  of 
the  energy  Indicated,  there  is  some  background  present. 
Furthermore,  the  low  energy  portion  is  hidden  due  to  back- 
ground pulses  (Compton  scattering,  bremsstrahlung,  random 
amplifier  and  detector  noise,  etc.).   To  determine  the  exact 
shape  of  the  spectrum  for  monoenergetic  gamma  rays  requires 
that  an  estimate  be  made  for  the  shape  of  the  low  energy 
region  (even  after  the  background  has  been  subtracted).   This 
estimate  is  needed  to  compute  the  total  number  of  counts 
in  the  spectrum. 

In  order  to  obtain  the  shape  of  the  spectrum  for  low 
energy  pulses  it  is  necessary  to  consider  the  different 
processes  that  have  been  employed  In  reproducing  a  spectrum, 
by  empirical  methods.   There  are  two  basic  procedures.   One 
is  the  Monte  Carlo  method  which  attempts  to  generate  the 
spectrum  by  tracking  individual  photons  through  a  crystal 
using  a  computer  (67,    68,  69).   Zerby  and  Koran  assumed  that 
the  source  was  In  a  vacuum,  thus  eliminating  background. 
Also,  the  calculations  are  first  carried  out  as  if  the  crys- 
tal produced  a  sharp  line  in  response  to  the  energy  deposited 
in  the  crystal.   The  response  spectrum  Is  subsequently 
spread  out  by  a  semi-empirical  Gaussian  broadening.   The 
physical  processes  that  gamma  rays  undergo  in  passing 
through  the  crystal  are  simulated  by  this  method.   The  Monte 
Carlo  calculation  Is  started  by  determining  the  direction 
of  incidence  of  the  gamma  ray  by  systematically  sampling 
the  isotropic,  monenergetic,  point  source  within  the  cone 
subtended  by  the  crystal.   The  different  interactions 
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(photoelectric,  pair,  Compton)  are  then  systematically 
sampled  In  relation  to  their  probabilities  of  occurring 

Zerby  and  Moran  argue  that  the  intercept  of  the 
spectrum  with  the  zero  pulse-height  value  can  be  calculated. 
They  note  that  photons  which. are  forward  scattered  only 
once  before  escaping  can  be  the  only  ones  that  contribute 
in  the  limit  of  zero  pulse-height.   They  obtain  an  expression 
for  this.   Unfortunately,  it  does  not  agree  well  with  the 
actual  zero-intercept  observed  in  carefully  performed 
experiments.   It  is  consistently  lower.   Thus,  some  difficulty 
in  determining  the  zero-intercept  exists. 

The  second  approach  is  to  employ  the  method  of  least- 
squares  to  fit  the  observed  spectrum  using  known  or  suspected 
energy  components  as  a  first  approximation  (70,  71).   These 
known  spectra  have  been  carefully  compiled  by  Heath  and 
others  using  stringent  laboratory  conditions  which  minimize 
the  effects  of  background  (72,  73).   The  idea  is  to  use  a 
program  to  fit  the  observed  spectra  using  a  complete  set  of  . 
functions.   It  is  then  desirable  to  be  able  to  adjust  the 
parameters  so  that  the  spectrum  response  for  any  energy 
gamma  ray  can  be  generated.   In  this  case  Heath  has  shown 
that  the  low  energy  region  can  be  approximated  by  a  linear 
extrapolation  to  the  zero  energy  axis  from  the  known, 
rather  flat  Compton  distribution.   For  the  present  experiment  a 
horizontal  extrapolation  has  been  employed. 
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